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Thia-paper  presents  a  detailed  exposition  of  the  problems  associated  with 
landing  high  performance  jet  fighter  aircraft  on  d  a  moving  aircraft  carrier 
along  with  a  new  procedure  for  designing  a  fully  h  matlc  flight  control  system 
tailored  to  the  aircraft  carrier  approach.  Indue  -n  the  development  is  a 
brief  description  of  the  present  SPN-42  ACI.S  (Aircraft  Carrier  Landing  System) 
and  the  probable  reason  for  its  very  limited  success.  Possible  alternatives 
to  the  SPN-42  system  have  suggested  a  more  sophisticated  ACLS  configuration,  in 
particular,  the  AMQAC  (Automatic  Multiloop  Optimal  Approach  Controller)  developed 
herein.  The  use  of  an  AMOAC  system  in  actual  carrier  approaches  should  yield  a 
substantial  improvement  to  present  automatic  carrier  landing  capability.// The 
results  of  Appendix  I  show  that  similar  closed  loop  dynamics  can  be  expec^d  for. 
other  carrier  based,  high  performance  aircraft.  In  particular,  it  is  felt  that 
the  AMQAC  design  procedure  has  potential  for  rapid  and  straightforward  application 
to  the  A-5A  and  the  F-4  aircraft. 
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PAKT  I 

INTRQDU JT ION 
A.  Background 

Aircraft  easier  landings  present  one  of  the  most  critical  problems 
currently  faced  by  Navy  pilots.  The  relatively  small  dimensions  of  the 
carrier  flight  deck  and  the  proximity  of  the  water  impose  severe  limi¬ 
tations  on  the  airplane's  allowable  position  errors.  The  motion  of  the 
aircraft  carrier  in  tne  sea,  pitch,  roll,  heave,  etc.,  further  complicates 
the  problem.  The  random  movement  of  the  carrier  flight  deck  and  the 
turbulence  created  by  the  carrier  interacting  with  the  free  air  stream 
both  demand  precise  control  of  the  airplane  during  the  carr:er  approach. 
High  performance  fighter/attack  aircraft  are  designed  primar? ly  for  op¬ 
eration  in  the  high  speed  portion  of  their  flight  envelopes;  the  lew  air 
speed  of  the  powered  approach  introduces  major  control  problems.  Ihe  ef¬ 
fectiveness  of  the  control  surfaces  is  severely  degraded  at  approach 
speeds.  Approach  air  speed  is  typicaiiy  below  minimum  drag  air  speed. 

The  reversal  in  the  slope  of  the  thrust  required  curve  introduces  right 
half  plane  zeros  into  the  sink  rate  to  control  surface  transfer  functions, 
thus  radically  different  control  laws  are  needed  in  the  approach  '-onfig- 
uration  than  in  the  remainder  of  the  flight  envelope.  The  need  for  im¬ 
proved  carrier  approach  control  is  borne  out  by  Naval  Air  Safety  Center 
statistics.  In  calendar  year  1965,  the  ratio  of  major  accidents  in  car- 
rior  approaches  per  total  flights  to  major  accidents  in  ground  approaches 
per  total  flights  was  approximately  ten  to  one.  In  calendar  v«ar  1967, 
this  figure  was  approximately  six  to  one. 


1 


NADC-SD-7153 

Past  efforts  in  applying  automatic  control  theory  to  the  carrier  ap¬ 
proach  problem  have  given  rise  to  the  Approach  Power  Compensator  (APC)  ana 
SPN-42  Automatic  Carrier  Landing  System  (ACLS).  The  APC  concept  is  based 
on  automatic  tnrottle  control  used  to  trim  thrust  during  a  power  approach. 

It  was  originally  designed  as  a  pilot  aid  for  manual  approaches.  For  this 
purpose,  it  has  achieved  some  success.  Subsequent  effort  has  been  given  to 
develop  an  automatic  elevator  controller  to  be  used  along  with  the  APC  to 
maintain  longitudinal  flight  path  control.  This  flight  path  controller  is 
designated  as  the  SPN-42  ACLS,  Mode  1.  The  SPN-42  system  has  been  under 
development  for  more  than  twelve  years  and  has  thus  far  achieved  very  lim¬ 
ited  success.  Difficulties  have  arisen  in  achieving  compatibility  of  the 
SPN-42  with  the  existing  APC.  The  SPN-42  is  presently  fleet  operational 
for  only  a  single  carrier  based  airplane  and  there  are  several  compatibility 
problems  with  other  Navy  airplanes. 

The  SPN-42/APC  combination  represents  a  reasonable  first  attempt  to 
design  an  automatic  carrier  landing  system.  The  carrier  approach  problem, 
however,  is  sufficiently  complex  that  it  requires  a  more  sophisticated  con¬ 
trol  system.  The  single  command  input  of  the  longitudinal  portion  of  SPN- 
42  elevator  deflection  and  the  thrust  trimming  of  the  APC  are  employed  to 
control  airspeed,  altitude  and  pitch  attitude  response  for  both  deterministic 
and  stochastic  inputs.  This  system  lacks  independence  in  specifying  the 
various  responses  for  both  inputs.  In  order  to  alleviate  this  problem  and 
obtain  good  response  characteristics  for  the  critical  variables,  a  multicon¬ 
troller,  multivariable  design  can  be  employed.  Reference  (a)  shows  the  ef¬ 
fect  of  multi loop  control  on  a  manual  carrier  approach. 

This  report  considers  a  longitudinal,  three  controller  system  incor¬ 
porating  conanand  inputs  of  elevator,  thrust  and  direct  lift  control  flaps 
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as  a  potential  advanced  ACLS.  The  control  system  configuration  is  a 
rather  general  12-parameter  feedback  system.  Design  procedure  is  based 
upon  the  Root  Square  Locus  technique  of  S.  S.  L.  Chang,  and  extended  to 
multi- input  systems  by  E.  Rynaski  and  R.  Whitbeck,  of  Cornell  Aeronauti¬ 
cal  Laboratory,  in  reference  (b).  For  the  purpose  of  this  report, good 
flying  qualities  will  be  defined  as  the  best  trade-off  between  (1)  mini¬ 
mum  dispersion  at  touchdown  and  (2)  smoothness  and  regularity  of  the  air¬ 
plane's  stability  and  control  characteristics,  and  riding  qualities.  In 
this  regard,  the  discussion  will  center  around  conventional  parameters  of 
the  characteristic  roots  of  a  dynamic  system. 


The  airplane  selected  for  this  study  is  the  F-8C  in  a  power  approach 
configuration.  Aerodynamic  parameters  and  physical  constants  for  this 
airplane  were  obtained  from  reference  (a). 

B.  Design  Considerations 

In  a  multivariable,  multicontroller  system  there  are  a  great  number 
of  possible  feedback  configurations.  The  determination  of  the  most  flex¬ 
ible  configuration  for  the  carrier  landing  problem  would  require  a  lengthy 
investigation  involving  at  least  12  independent  feedback  parameters;  that 


is,  a  feedback  loop  from  each  measured  output  6,  d,  u,  and  h  to  each  con¬ 


trol  input  6g,  6p,  and  <5^.  Such  an  investigation  would  be  extremely  com¬ 
plicated  and  tax  the  capabilities  of  current  computing  equipment.  For 
this  reason,  a  method  was  sought  to  evaluate  a  12-parameter  feedback  con¬ 


figuration  as  a  function  of  a  few  independent  parameters;  namely,  four 
independent  variables. 

To  this  date,  research  on  optimal  control  systems  has  shown  that 
for  linear  multidi'mensional  problems  the  optimal  feedback  configuration 
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tends  to  produce  closed  loop  systems  which  compare  favorably  with  conven¬ 
tional  design  criteria.  With  the  introduction  of  optimal  control  tech¬ 
niques,  the  number  of  independent  parameters  has  been  reduced  to  those 
required  to  specify  the  performance  index. 

Optimal  control  techniques  dictate  the  form  of  the  performance  in¬ 
dex;  namely,  quadratic  forms  in  the  controls  and  the  measured  outputs. 

PI  =  1/2  J  [  a;  u2  *  a2  h2  +  a3  02  +  a4  02 

+  di  62;2  +  d?  6f2  +  d3  6T2  -I  dt 

In  adopting  this  performance  measure  to  the  carrier  landing  pro¬ 
blem,  some  simplifications  based  on  •.  sical  control  characteristics  can 

be  made.  The  shallow  descent  angle,  which  requires  more  precise  control 
•  • 
of  h  than  u,  allows  aj  =  0  while  the  redundance  of  controlling  0  and  0 

allows  a3  »  0. 

Upon  normalizing  by  d?,  this  measure  becomes- 

Joo  ^  ^  ^ 

[  k]  h2  +  k2  e2  +  d[  6k2+  d2  6p2  “  dt 

The  optimal  feedback  control  law,  that  is,  that  feedback  con¬ 
figuration  which  minimizes  PI  can  now  be  evaluated  as  a  function  of  the 
four  parameters  ki,  k2,  dj  and  dz.  The  remaining  problem  is  to  characterize 
the  flying  qualities  of  the  resulting  closed  loop  air  frame  equations  as 
the  four  parameters  of  PI  are  varied.  This  is  accomplished  via  the  Root 
Square  Locus  technique  of  reference  (b). 
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PART  II 


THEORY 

A.  Optimal  Control  Equations 

The  optimal  system  is  described  by  three  differential  equations. 

The  first  equation  is  the  constraint  imposed  by  the  basic  air  frame  longi¬ 
tudinal  equations  of  motion.  The  second  and  third  equations  are  the 
Euler  variational  equations  of  the  optimal  system  which  are  obtained  by 
using  the  calculus  of  variations.  A  fourth  equation,  a  matrix  Ricotti 
equation,  determines  the  feedback  gains. 

Air  Frame  Longitudinal  Equations  of  Motion  (Stability  Axes) 

Figure  1  defines  the  stability  axes  system 

Drag:  u  =  X,  u  ♦  X  w  -  (gcosyo)0  +  X,.,.  6_ 

vJ  ^  Oc  L 

+  X6T  +  X^p.  6p  (la) 

Lift  w  »  Zy  u  +  Zw  w  +  Vo  6  (g  sin  y0)6 

+Z6E«E  +  ’Z6T  «T  +  Z6F  6F 

'  *  f 

Pitching  Moment  0  *  M  u+M  w+M'w+Mq  6+  M.p  6 

u  w  w  oc  t 

+  j«6F  <5p  ♦  M6t  6t  (lc) 


Altitude 


Vo  e  =  w  +  h 


where  u,  w,  h,  0,  0,  6^,  6p  and  6^  are  incremental  quantities  and  the 
aerodynamic  parameters  are  as  defined  in  Appendix  A. 

For  the  initial  portion  of  this  report,  actuator  lags  and  non-  • 
linearities  will  be  neglected.  In  later  sections,  these  effects  will  be 
introduced. 

Equations  (la,  b,  c,  d)  can  be  rewritten  in  matrix  notation  as: 
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*Xw  XwV0^cosY0 

Zw  gsinYo  -  2WV0 

0  0 


-(M  ♦M'.Z  Vo  +  M'Z  Vo  * 
v  w  w  vr  w  w  w 


-gMwsinY0 


A  more  compact  form  of  equation  (2)  is 

x  =  Ax  +  Bj5  (3) 

•  » 

where  <  =  column  (u,  h,  8,  9),  6  -  column  (.6  , e>  '  and  A  and  B  are  the 

r  I  r 

corresponding  matrices  of  equation  (2). 

After  substituting  the  values  of  aerodynamic  constants  from  Appendix  A, 

A  and  B  become: 
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- 

-0-0506 

0.0073 

-33  8 

0 

0  2o4 

-0.554 

127  4 

0 

0 

0 

0 

1  0 

-0  000583 

0.00465 

-1  051 

-0.390 

-0.672 

0  0001595 

-6.72 

19  95 

0  0002104 

i8 . 9 

0 

0 

0 

2  196 

-4  67  X  10'6 

0  0357 

Euler  Eqctions 

The  carrier  landing  problem  can  be  formulated  as  follows: 

Given  the  air  frame  constraint  equation  (3),  find  S  as  a  linen,-  function 
of  x_,  linear  feedback,  such  that  the  Pi  is  minimized  for  values  of  k', 
k2,  bi.  and  b?  which  result  in  good  airplane  flying  qualities 


Expressing  PI  in  matrix  notation  yields: 

PI  *  1/ 2^  [x  rK  x  +  6  ru  ±  1  dt 


0  0  0  0 
a 

where  K  =■  0  kj  0  0 


di  0  0 

where  K  *  0  kj  0  0  and  D  =  010 

0  0  0  0  0  0  d2 

0  0  0  k' 

x 

-  d  L  J 

Since  the  solution  of  this  problem  is  constrained  to  satisfy  equa¬ 
tion  (3),  a  vector  function  of  time  £(t),  the  aajoint  variable  is  intro¬ 
duced.  Introduction  of  the  adjoint  variables  allows  the  constraint  equa¬ 
tion  to  be  included  under  the  integral  in  equation  (4)  without  changing 
the  value  of  PI. 

PI  *  1/iJ^  [x  rK  *  *  6rD6  ♦  £  r  (-X  »Ax  ♦  B6JJ  dt  (S) 
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but  'i.  +  Ax  +  B6  =  0 

or  PI  *  PI 

Using  the  results  of  Appendix  B  yields  equations  (6a  and  f)» 

£  ■=  -A  T£  -  Kx  (6.  ) 

t)^  =  -B  T  £  (6b) 

Equations  (6a  and  b)  are  the  Euler  equations  of  the  optimal  system. 
Equations  (6a  and  b)  along  with  equation  (3)  comprise  the  necessary  con¬ 
ditions  for  an  optimal  system  Using  a  matrix  format  and  eliminating 
equation  (6b)  yields1 


Feedback  Gains 

Once  a  set  of  parameter  values  is  determined,  it  is  then  necessary 

i 

to  determine  the  corresponding  feedback  gains  Equation  (6b)  indicates 
that  the  optimal  control  £(t)  is  a  linear  combination  of  the  adjoint 
variables  £(t) 

6_(t)  =  -U':b‘  £  (t)  (8) 

In  order  to  obtain  -S(t)  as  a  linear  combination  of  x_(t),  or  linear 
feedback,  it  is  required  that  the  foil  wing  relationship  be  satisfied: 

£(t)  =  P(t)  x(t)  (9) 

In  which  case: 

6(t)  =  -D-1BT  P(t)  x (t )  (10) 
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Appendix  B  shows  that  equation  (9)  is  indeed  valid  and  also  that 
the  infinite  upper  limit  of  integration  in  PJ  guarantees  that  P(t)  is  a 
ma'ri*  of  constant  elements  In  Appendix  C  an  algebraic  equation  for  P 
is  derived 


0  =  PA  *  ATP  *  K  -  PBD  BfP  (11) 

In  this  investigation,  the  solution  of  equation  (11)  was  accom¬ 
plished  with  a  digital  computer  program  developed  by  Cornell  Aeronautical 
Laboratories  This  program  implements  a  scheme  based  upon  the  eigen 
vector  expansion  of  the  Hamiltonian  matrix  This  method  was  developed 
by  Potter  and  was  first  published  in  reference  (c)  The  P  matrix  is  then 
substituted  into  equation  (10)  to  obtain  the  feedback  gain  matrix,  H 

After  computation  of  the  feedback  gams,  the  overall  closed  loop 
system  can  be  examined  This  is  accomplished  by  calculating  the  transfer 
functions  of  output  variable,  to  the  command  inputs  Starting  with  the 
open  loop  air  frame,  equation  (3),  and  the  equation  of  the  feedback  con¬ 
troller  which  is 

i  '  £c  ♦  if  b  02) 

where  Ac  is  the  command  input  and  5^  .  is  the  feedback  controller 

output  given  by 


the  closed  loop  equations  become: 

x  -  (A  -  BH)x  +  B6 c 

Laplace  transforming  equation  (14)  and  solving  for  X(s)  yields 

X(s)  *  (Is  -  A  ♦  BH]'  B  A(s) 


(13) 


(14) 


(15) 
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Pro*  equ*non  (15)  it  is  found  tint  tin-  closed  loop  txansfer  matrix 

Mc'.(s)  --  f  Is  A  -  BH)"  B  (lb) 

B  iLcaim.  Procedure 

fn  order  to  determine  which  values  of  k,  ,  ,  d-  and  d?  result 

in  good  flying  qualities,  the  closed  loop  roots  of  the  optimal  system  are 
determined  as  a  function  of  these  parameters  via  a  digital  program  im¬ 
plementing  th*1  Hoot  Square  Locus  technique  The  program  determines 
a  region  ol  parameter  values  with  good  short  period  and  phugoid  charac- 
t or i st ! i s  Phis  process  is  discussed  n  detail  in  section  II  B-l  and  in 
Appendix  (Cl  The  resulting  sreadv  state,  matrix  Kuotti  equation  is  then 
solved  and  the  complete  set  of  multi  loop  feedback  gains  are  evaluated 
In  this  manner,  the  use  of  optimal  control  theory  has  transformed  the 
problem  of  choosing  twelve  feedback  gains  to  a  lesser  problem  of  choosing 
four  performance  index  parameters  The  design  procedure  then  proceeds 
along  a  more  conventional,  state  variable  format. 

Since  the  Hoot  Square  Locus  has  been  employed  only  to  achieve  satis 
factory  closed  loop  system  roots,  it  is  necessary  to  investigate  the  re¬ 
sulting  closed  loop  numerator  zeros  Thus,  the  closed  loop  transfer  matrix 
is  calculated  and  its  component  terms  invent ig^'ed  along  the  guidelines 
established  in  section  li  C- l  If  the  numerator  terms  cannot  mef  these 
guidelines,  the  performance  index  parameters  are  adjusted  either  by  re- 
employing  tne  Root  Square  Locus,  or  on  the  basis  of  more  immediate  intel¬ 
ligence.  Upon  reaching  a  satisfactory  closed  loop  transfer  roatr,v,  those 
preliminary  feed  back  gains  would  be  incorporated  into  a  carrier  approach 


simulation. 
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This  simulation  is  exercised  to  evaluate  the  operation  of  the  AMOAC 
control  system  in  an  accurate  model  of  the  carrier  approach  environment, 

4  |  '  , 

This  is  accomplished  as  follows: 

a  The  system  is  first  operated  independently  of  any  glide  , 
slope  positional  error  information,  as  an  inner  loop  augmented  aircraft, 

I  i 

i 

without  actuators  or  aircraft ,non-linearit ies  and  lags.  This  is  done  to 

■  .  '  '  i 

establish  the  aircraft's  response  to  airstream  velocity  disturbances  or 
steady  state  aircraft  carrier  burble  and  .pitch,  induced  burble 

I 

b  The  aircraft's  control  surface  and  propulsive  non-linear¬ 
ities  and  lags  are  incorporated  into  the  tlosed  loop  system  and  the  re¬ 
sulting  simulation  exercised  again  as  an  inner  i  loop  augmented  aircraft. 

l 

c  A  glide  slope  positional  error  loop  of  the  form 

hffe  *  h  +  k  h  f  1 7) 

or  upon  Laplace  transforming  \ 

Hf b (* )  »  fl  ♦  £  jH  fs)  '  (  (18) 

is  closed  around  the  augmented  aircraft  of  paragraph  b  (above'i.  Exer¬ 
cising  this  simulation  demonstrates  the  capacity  of  the  AMOAC  control 
system-  to  effect  accurate  glide  slope  control  when  combined  with  the 
simple  outer  navation  loop  described  by  equation  fl8J 

d.  The  closed  loop  system  of  c  is  then  tested  along  a  turbu- 

I 

lence  free  glide  slope  toward  a  desired  TUP  (touch  down  point'  whose  co¬ 
ordinates  vary  randomly  in  a  manner  appropriate  to  an  aircraft  carrier 
in  severe  sea  swell  conditions  The  resulting  vertical  touch  .down  dis¬ 
persion  fo^)  is  calculated  from  the  carrier  motion  power  spectral  densi¬ 
ties  and  the  closed  loop  transfer  matrix 
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If  an  unsatisfactory  result  is  obtained  from  any  of  the  four  simula¬ 
tion  exercises,  the  performance  index  parameters  are  adjusted  and  the 
design  procedure  re-employed.  After  satisfactory  results  are  obtained 
from  simulation  d,  the  AMOAC  system  design  is  considered  satisfactory 
and  the  design  procedure  is  terminated.  Figure  2  demonstrates  a  func¬ 
tional  block  diagram  for  this  design  p  ocedure. 

Figure  2  breaks  down  the  design  process  into  two  phases-  The  first 
phase,  i,e.,  system  synthesis,  employs  optimal  control  theory  to  obtain 
a  preliminary  design  configuration.  Conventional  state  variable,  analysis 
techniques  are  employed  in  the  second  phase  in  order  to  evaluate  the  pre¬ 
liminary  design. 


This  second  phase  is  broken  up  into  five  separate  analyses,  each 
of  which  must  be  satisfied  before  the  design  is  accepted.  After  each  of 
the  five  evaluation  steps,  a  decision  is  made  whether  to  continue  to  the 
succeeding  step  or  adjust  the  multiloop  gains.  The  evaluation  procedure 
is  broken  into  relatively  small  increments  in  the  hope  that  only  small 
adjustments  of  the  performance  index  parameters  will  be  required  to  re¬ 
tune  a  preliminary  design  configuration.  The  results  of  this  investigation 
appear  to  validate  this  design  approach, 

1 .  System  Koots 

The  root  locations  of  the  closed  loop  optimal  system  vary  as 
functions  of  the  performance  index  parameters.  This  section  describes 
how  these  root  locations  are  determined.  The  root  square  locus  technique 
described  in  Appendix  (C)  is  now  adapted  for  use  in  the  AMOAC  design  pro¬ 
cess.  In  Appendix  C,  the  root  square  locus  equation  was  derived  and  is 
slightly  modified  here  for  convenience. 


DESIGN 

EVALUATION 


CONVENTIONAL 

STATE 

VARIABLE 

ANALYSIS 


FIGURE  2,  AMOAC  LOGIC  FLOW 
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det  (I  -  D-;MT  (-$)  K^M  ( s )  )  -  0  (19) 

The  elements  of  MCs)  are  designated  as  M.  (s) ,  where  the  first 

j 

index  denotes  the  measured  output,  and  the  second  index  the  control  input. 

t 

Sub*-t_ »-ut ing  for  D  and  K  in  equation  (19)  yields: 

A  A  ■  '  h  ---A 

0  =  1  ♦  k-  M;>/  +  k?  M.?  •  k(  bj  ♦  k?  bi  M,,j  +  k,  k?  bA  M^2  ♦ 

*  *  A  .  .  A  *  *  A  *  •  •  ^ 

♦  k  b  M  ♦  k  b  M  '  U  1)  M'.i  *  k,  k.  b  b,  M*s2  f?0) 

where  A.  -  (s)  (-si,  M?.  -  Mj  M?j  and  bj  =  1/d*  . 


•Hic.h  element,  Mi-(sJ.  of  M(s)  is  of  the  form 

N.4(s) 


Mij(S)  =  Sts 


where  ;  (s)  is  the  characteristic  polynomial  of  the  open  loop  system,  lhe 

<•?  ? 
terms  of  tlie  form  Mj^  appear  to  have  denominators  of  (£(s))  •  however, 

reference  (b)  has  shown  that  the  numerator  of  such  terms  includes  a  factor 


of  A(s). 


N  42 

l  .  A (s  i 


(.  d  (s)2) 


Rearranging  equation  (20)  yields  the  expression  used  in  this  report 


to  determine  the  closed  loop  roots. 


Q  s  L  +  k>  (•  (ft?;  *  b:  ft?i  ♦  b2^2  3)  1  k2 

&(s)  *  k?  (ftc;  *  bj  Ne, ;  +  b2  ftu?) 


'  *  ■*  7  '  ■'  '  A.  , 

fl> l  fti2  *  b?  jki  ♦  b:bi  jjil)) 


15 


NADC-SD-7153 


Equation  (21)  has  form  of  a  conventional  multi  loop,  root  locus.  It 
consists  of  an  outer  loop,  closed  around  parameter  k.,  a  numerator  and 
denominator  loop  each  closed  simultaneously  abound  k?  and  three  inner 
loops  with  b_  and  b2  as  parameters. 

The  three  inner  loops  differ  slightly  with  the  conventional  root 
locus  format.  They  are  each  functions  of  the  same  two  parameters  bj 
and  b;  Because  of  the  b;  b 2  term,  they  cannot  be  broken  down  into  two 
single  parameter  loops.  Each  of  the  inner  loops  has  a  distinct  physical 

interpretation.  The  first  inner  loop  of  the  numerator  is: 

A  -A  .A 

0  «  N22(s)  ♦  b|N?,  (s)  ~  b?  K23fs)  (22) 

and  contains  the  numerators  of  all  the  sink  rate  transfer  functions.  This 
loop  weights  the  importance  of  elevator,  thrust  and  flaps  as  primary  con¬ 
trollers  of  sink  rate.  This  loop  will  be  referred  to  as  the  sink  rate 
control  loop. 

The  inner  loop  of  the  denominator  is: 

A  .  A  -A 

0  ^  N^Cs)  +  bi  Ni,  x  (s)  *  b?  N„3  (s)  (23) 

and  consists  of  the  numerators  of  the  pitch  rate  transfer  functions  and 
is  thus  called  the  pitch  rate  control  loop  This  loop  weights  the  con-  . 
tribution  of  the  three  control  inputs  as  primary  controllers  of  pitch 
rate. 

The  second  inner  loop  of  the  numerator  is: 

.A  ,  .  .  A 

0  =  bt  NiJ  (s)  +  b2  (s)  ♦  b.  b;.  N^?(s).  (24) 

and  contains  the  cross  coupling  terms.  This  loop  takes  into  considera¬ 
tion  the  interaction  between  any  two  of  the  three  control  actions  and 
is  referred  to  as  the  cross  coupling  loop. 
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Equation  (21)  is  an  eighth  order  polynomial  in  s,  with  roots  that 
are  symmetric  about  the  ju>  axis.  Each  term  m  equation  (21)  is  the  pro¬ 
duct  of  a  polynomial  P(s)  and  another  polynomial  P(-s).  The  roots  of  eq¬ 
uation  f 2 1 )  correspond  to  the  closed  loop  roots  of  equation  (7),  i,e  . ,  the 
roots  of  the  optimal  system  and  the  adjoint  system.  If  values  of  the 
parameters  are  restricted  to  be  positive,  the  R  matrix  of  equation  (4)  is 
positive  definite  and  the  Q  matrix  positive  semidefinite  Under  thi* 

restriction,  PI  becomes  a  Liapunov  function  for  the  closed  loop  system  and 

*  »■  •  » 

those  open  loop  roots,  k  =k  =b  =  b  =0,  in  the  left  half  plane  must 

I  ?  1  / 

remain  in  that  plane  as  the  parameters  are  increased 

Since  the  basic  airframe  has  stable  roots,  the  part  of  the  root 
locus  in  the  left  half  plane  must  correspond  to  the  closed  loop  optimal 
system 

The  design  procedure  followed  in  this  investigation  was  to  first 
*  > 

determine  values  of  b;  and  b2  that  provide  roots  for  shaping  good  numera¬ 
tor  and  denominator  loops..  Then  ki  and  k?  was  chosen  to  produce  satis¬ 
factory  phugcid  and  short  period  roots 

In  order  to  obtain  the  roots  of  equations  (22)  through  124),  a 
digital  program  was  set  up  to  find  roots  of  an  equation  of  the  following 
form. 

0  =  (J)  ♦  C1  (s)  +  e2  (s)  ♦  C;  c2  K  (s) 

*  *  * 

This  program  is  then  used  to  obtain  proper  values  of  kt,  k2,  b  and 
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?.  Carrier  Approach  Simulation 

The  heart  of  the  AMOAC  design  evaluation  is  the  exercise  of 
simulated  carrier  approaches  via  digital  computer  programs  which  accu¬ 
rately  model  the  carrier  approach  environment  and  the  closed  loop  airplane 
dynamics  This  section  describes  the  carrier  approach  environment,  the 
environmental  model  used  in  the  simulation,  the  simulation  itself,  and 
its  exercise  ground  rules, 

a.  Approach  hnvironment 

The  typical  problem  of  flight  mechanics  is  to  determine  the 
motion  of  an  aircraft  with  respect  to  its  surrounding  airstream,  when  it 
is  subjected  to  some  excitation  forces.  The  longitudinal  equations  o. 
motion,  equation  (3J ,  is  derived  in  Part  II,  Section  A  m  accordance  with 
this  format.  When  this  procedure  is  used  to  investigate  a  conventional 
powered  landing  approach,  however,  a  new  problem  arises.  This  problem  is 
the  referral  of  the  aircraft's  coordinates  to  the  inertial  coordinates  of 
the  desired  TDP,  (touchdown  point)-  If,  in  addition,  the  TDP  is  not  sta¬ 
tionary  with  time,  the  problem  is  further  complicated,  A  final  complexity 
is  the  non-uniformity  of  the  airstream  velocities.  All  of  these  factors 
together  describe  the  carrier  approach  problem. 

The  first  consideration  in  describing  the  carrier  approach 

environment  is  the  undisturbed  glide  slope.  When  the  airplane  is  in  the 

final  approach,  it  tries  to  fly  a  straight  line  flight  path  to  the  TDP- 

In  this  process,  a  constant  airspeed,  VOJ  and  a  constant  angle  of  attack, 

a„.  .  are  maintained.  This  airspeed  is  essentially  the  minimum  safe  air-- 
1  run  1  7 

speed  for  the  airplane  based  upon  stall  and  other  aerodynamic  considerations, 
and  is  determined  by  V0  and  the  airplane’s  weight,  W  The  value  of 
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sink  rate,  h0,  and  consequently  the  glide  slope  angle,  Yo, 


Y  o  85  tan"  1  (ho/V0  ) 


(25) 


are  set  by  structure,  and  safety  limits  at  deck  contact.  The  values  for 
these  parameters  were  obtained  from  reference  (a)  and  are  listed  in  Ap¬ 
pendix  A. 

The  aircraft  carrier,  however,  maintains  a  velocity  in  the 
same  positive  x  direction  as  the  airplane.  An  equivalent  situation  is  a 
stationary  carrier  with  the  airstream  moving  past  the  carrier  and  parallel 

to  it.  In  this  situation,  the  closure  rate  of  the  airplane  and  the  TDP, 

* 

Vo,  is  (Vo  -  WOD),  where  WOD  stands  for  wind  over  deck,  namely,  aircraft 
carrier  forward  speed.  Figure  3  shows  this  situation.  Note  that  angle 
of  attack  and  flight  path  angle  are  now  calculated  with  respect  to  the 
new  closure  rate.  t 

Y  o  -  tan'1  (ho/V0  )  *  -4  deg.  (26) 


aTrim  *  aTrim 


*  (Y°  -  Yo)  =  7.4  deg 


where  VQ  equals  180  feet  per  second.  The  reference  system  of  figure  3 
with  its  origin  at  the  TDP  of  the  stationary  carrier,  will  be  used  through 


out. 


As  the  airstream  passes  over  the  carrier,  the  flow  is  severely 
disturbed.  This  gives  rise  to  a  flow  pattern  which  varies  with  the  X 
and  Z  coordinates  of  figure  3  but  is  stationary  with  time.  The  influence 
of  sea  waves  imparts  a  random  movement  to  the  TDP  which  further 
disturbs  the  airstream.  Figure  4  illustrates  the  carrier  approach 
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envirortment . 

The  model  of  the  approach  environment  was  broken  into  two  parts. 
The  first  part  is  a  model  of  the  airstream  variation  for  points  along 
the  undisturbed  glide  slope.  The  resulting  geometry  is  shown  in  figure 
S.  Carrier  burble,  i.e.,  the  perturbed  portion  of  the  airstream  is  de¬ 
composed  into  a  horizontal  component,  Ug,  and  a  vertical  component,  W 
each  of  which  is  a  function  of  the  range,  R.  The  results  of  referenced 
and  e  indicate  the  following  form  for  carrier  burble: 


YR> 


-2d.0«  ,  17. Se  . R  <  2C00  ft. 

*23. Oe  ,  17.3e-R/75  a 

+  (1.8  -0.G009R)Sin(26.0  -  G.024R).R  ->2000  ft 


YR) 


cos  (R/255) . , ,  r  <2306  ft. 


C  7  naa  -R/11S4 
I  ■7-0e  COS  (R/255).,. 

■  -7.0e  R/1154  cos  fR/255j  + 
V^+(4.15-  0-0185R)Sin(26.0  -  0.024 


018SR)Sin(26.0  -  0.024R)R  <2306  ft 


This  model  is  quite  conservative.  The  maximum  values  of  u  .and 

\  are  14  feet  per  second  and  11  feet  per  second  respectively.  Both 

maxima,  values  occur  nearly  midway  between  the  ramp  and  the  TOP. 

The  second  phase  of  the  approach  environment  model  consists  of 

a  description  of  the  random  carrier  motion.  Reference  d  indicates  the 

for.  of  typical  power  spectral  densities  for  an  aircraft  carrier  in  pitch, 

ecm<  hMV'<  2cm  ,nd  sur*e’  *cm  for  *  severe  sea  state.  The  particular 
model  used  in  this  report  is: 


FIGURE  5.  CARRIER  APPROACH  MODEL  (STATIONARY  TD) 


0.0  feet 


(30a) 
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Surge ,  x  x 

6  cm  cm 


=  £  [x  J 
1  cnr 


=  £  [  (x  -  x  )'■  \ 
1  cm  cm  ’ 


=  1.0  feet7  (30b) 


0.2023fajj _ _ 

0.05883a)2  +  (0.5776  -u)2)7 


(ode ) 


Heave,  z  z  m  =  £  (z  ] 
cm  cm  v  cmJ 


*  0.0  feet 


(31a) 


o z7  =  H  I ( 2  -  z  )7J=2.0  feet" 

cm  cm  cm  J 


0. 4198a)7 
0.01795a)2 


Pitch,  0  0  =  £  [0  1  =  0.0  degree 

cm  xm  cm'  ° 


(31b) 


♦  (0.5476a)7)2  (31c) 
(32a) 


=  £  re  -e  y  J 
1  cm  cm'  J 


•±  1.0  degrees  (32b  j 


“ci-)  * 


0.2217a)7 

0.06S16a)7 


♦  (0.6084  -  a)2 ) : 


(32c) 


Where  Sx (u>)  denotes  the  power  spectral  density  of  x  and  is  defined  as: 


..  ,  ,  I  -ia)T,lim 

V*> s  1  c 


^fT 

;,T  :  -T 


X(t)x(t-T)  dtj  dl 


b.  Approach  Simulator 

The  approach  simulation  is  also  a  two  part  procedure.  The  first, 
part  evaluates  the  response  of  the  augmented  aircraft  to  the  carrier  dis¬ 
turbed  airstrcam  as  modelled  in  equations  (28)  and  (29).  Terminal  errors 
due  to  carrier  deck  motion  art  determined  in  the  second  phase  of  the  ap¬ 
proach  simulation. 


:  •>  JTV  ^ 
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i.  Carrier  Burble  Response 

Figure  5  is  an  accurate  representation  of  the  geometry 
and  environment  model  used  in  this  digital  simulation.  Tho  axis  system 
originates  at  the  stationary  TP  point  with  the  Z  axis  in  the  positive 
vertical  direction  and  X  axis  along  the  flight  deck  pointed  aft.  The 
nominal  glide  slope  emanates  from  the  TDP  inclined  by  the  angle  v*  - 

O 

4  degrees  from  the  X  axis.  The  carrier  burble  along  the  glide  slope  is 
characterized  by  two  functions  of  range,  R,  according  to  equations  (28) 
and  (.29). 

The  aircraft’s  linearized  equations  of  motion  are  solved 
about  the  nominal  glide  slope  to  yield  instantaneous  values  of  sink  rate 
and  airspeed,  relative  to  a  stationary  airstream.  These  values  are  then 
integrated  and  added  to  the  nominal  glide  slope  to  yield  the  aircraft's 
position.  The  carrier  burble  is  then  computed  and  incorporated  into 
sink  rate  and  airspeed  and  these  quantities  used  in  evaluating  the  equa¬ 
tions  of  motion. 

The  perpendicular  distance  between  the  nominal  glide  slope 
and  the  aircraft's  position,  h£  is  then  generated.  This  position  error 
is  multiplied  by  a  preselected  gain  and  added  to  the  sink  rate  feedback 
term.  This  procedure  simulates  a  crude  navigation  command  based  upon 
shipboard  telemetry.  This  navigation  guidance  is  most  likely  not  the 
best  available  but  it  can  demonstrate  the  performance  capability  of  AMOAC 
in  a  closed  loop  ACLS  (aircraft  carrier  landing  system) 

Figure  6  illustrates  the  geometry  used  in  generating  the 
aircraft's  flight  path.  An  accurate  description  of  the  simulation  me¬ 
chanics  is  given  in  Appendix  D. 
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-FIGURE  6,  CARRIER  APPROACH  SIMULATION  GEOMETRY 
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One  important  aspect  of  the  aircraft's  response  is  its  verti¬ 
cal  clearance  of  the  aft  extreme  of  the  flight  deck,  namely,  its  height 
over  the  ramp.  HR.  second  in  importance  only  to  the  TOP  dispersion,  AX,,. 
However,  due  to  the  point  mass  nature  of  the  translational  dynamics  and 
the  lack  of  hook  and  cable  contact  geometry  in  this  simulation,  the 
considerations  were  replaced  by  a  vertical  error  tolerance  over  the  TOP, 
«V  In  this  manner,  it  is  possible  to  assess  the  flight  path  response 
of  AMOAC  in  a  realistic  environment,  with  a  minimum  of  detailed  geometry, 
in  evaluating  the  performance  of  the  augmented  aircraft,  an  unsatisfactory 
error  in  HR,  41^,  was  chosen  to  be  any  excess  of  one-half  of  the  nominal 
H„.  The  value  of  »HU  was  chosen  to  be  the  vertical  distance  along  the 
glide  slope  that  corresponds  to  the  horizontal  distance  between  two  co„- 
aecutive  arresting  cables.  For  modern  carriers  XR  is  approximately  200 

feet  with  40  feet  between  arresting  cables.  Thus,  satisfactory  performance 
requires : 

rHJ  -  7-° fget  (34) 

|aHu|  i  1.8  feet  (3S) 

ii.  Carrier  Motion  induced  Dispersion 

This  computation  yields  the  mean  square  error  in  yerti- 

cal  position  over  the  TOP  induced  by  random  motion  of  the  aircraft  car: 

rier  according  to  equations  (30),  (31)  and  (32).  The  necessary  equations 

in  Appendix  £  assuming  shipboard  telemetry  provides  the  ship's 

motion  infomation  to  the  aircraft  control  system.  The  following  is  a 

condensation  of  the  results  of  Appendix  £  It  ic  nn 

c*  ls  rirst  necessary  to  intro¬ 
duce  an  augmented  state  variable  * 


. . ^i,»  riJI 
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£  s  column  (x,  h.  u,  h,  0,0  ) 

Associated  with  x  is  the  foliowin,  equation  of  motion 


where 


'V  -V  -v,  n 

X  *  A  x  ♦  B  6 


0  0 


0  0 


0  0 


0  0 


0  0 


0  0 


ooo 

ooo 


i  ooo 

0  10  0 


and  the  control  law 


where 


H  *  [HJ 


0  0  10  0 
o  K  o  1  o 
0  0  0  0  1 


If  °h  denotes  the  mean  square  error  in  vertical  aircraft  position 


due  to  carrier  motion,  then 
where  c  h2  *  j  s.  („ 


*  S»W<i“  (39, 

\M  .  |A22Mp  s2f  M  -  |A25Mt  |s 

m  wc 

AtuO  -  [jui  -  &  ♦  8HJ* 1  [juiI-Aj 


J 


■’".MIT 


IS!* 


!W 
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A  digital  program  performs  the  necessary  computations, 
c.  Simulation  Exercises 

In  the  evaluation  of  a  preliminary  AMOAC  design,  the  car¬ 
rier  approach  simulation  is  exercised  several  times. 

At  first,  a  simple  aircraft  model  is  used,  i.e.  no  control 
actuators,  lags  or  nonlinearities,  and  there  is  no  navigation  loop  clo¬ 
sure,  i.e.,  no  hg  feedback.  This  run  is  used  to  determine  the  dynamic 
response  of  the  aircraft’s  variables  and  controllers  to  the  carrier  ap¬ 
proach  environment.  The  ability  to  control  the  aircraft’s  rate  responses, 
i.e.,  sink  rate,  airspeed,  flight  path  and  pitch  rate  as  well  as  the  var¬ 
ious  control  expenditures  are  the  prime  concern  of  this  run.  If  it  is 
discovered  that  a  variable  is  inadequately  controlled  or  a  controller  is 
used  too  much  or  too  little,  its  respective  performance  index  parameter 
can  be  adjusted.  The  root  square  locus  is  then  used  to  determine  a  new 
preliminary  design. 

The  next  exercise  of  the  approach  simulation  evaluates  the 
effects  of  controller  actuators,  lags  and  non-linearities.  The  aircraft 
model  is  updated  to  include  these  controller  dynamics  and  non-linearities, 
and  the  carrier  approach  is  run  and  evaluated  along  the  same  ground  rules 
as  the  previous  run. 

After  a  preliminary  design  is  checked  out  successfully  in 
the  previous  two  exercises,  the  AMOAC  control  system  is  incorporated  into 
the  ACLS  system  described  in  the  previous  section  and  the  simulation  is 
again  exercised.  In  this  case,  two  runs  are  made.  The  first  run  utilizes 
a  twenty  foot  initial  vertical  glide  slope  error  and  no  burble.  A  success¬ 
ful  design  is  called  upon  to  null  out  all  errors,  i.e.,  position,  rates 


I 
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and  controller  deflections,  within  10  seconds  or  1800  feet  along  the  glide 
slope.  For  the  purpose  of  this  simulation,  nulling  out  errors  refers  to 
reducing  initial  errors  to  10  percent  and  maintaining  that  level  or  less. 
Airspeed  and  sink  rate,  which  are  initially  zero,  must  be  maintained  with¬ 
in  one  foot  per  second  and  one-half  foot  per  second  respectively  for  an 
error  to  be  called  null. 

The  aircraft  is  then  started  20  feet  off  the  nominal  glide 
slope  and  the  closed  loop  ACLS  is  allowed  to  land  the  aircraft  in  the 
presence  of  the  burble.  The  AH^  and  AH^  requirements  are  then  checked. 
Successful  runs  lead  the  computation  of  carrier  motion  induced  dispersion 
which  completes  the  simulation  exercises. 

The  choice  of  a  guidance  loop,  i.e.,  closed  loop  ACLS,  is 
accomplished  with  the  aid  of  the  analysis  of  the  next  section. 

3.  Guidance  Loop  Closure 

In  the  process  of  evaluating  the  AM0AC  system's  ACLS  capabili¬ 
ties,  the  introduction  of  a  guidance  loop  became  necessary.  The  guidance 
scheme  was  chosen  to  be  a  simple,  single  position  feedback  loop.  In' par¬ 
ticular,  the  perpendicular  deviation  from  the  nominal  glide  slope,  i.e., 
h.  of  figure  6  was  woighted  with  a  gain  k,  and  added  to  the  h  feedback 


hfb  *  h  *  k  he  « 

For  the  purposes  of  analysis  and  selection  of  a  proper  value 

for  k,  the  hg  term  was  considered  approximately  equal  to  the  vertical 

glide  slope  error.  For  a  4  degree  slope,  the  resulting  error  is  less 

than  one-half  of  one  percent.  With  this  approximation,  equation  (40) 
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becomes : 


hfb  •  h  *  k 


In  appendix  F  equation  (41)  is  incorporated  into  the  optimal  closed  loop 
system  equations,  equation  (3)  and  13),  and  the  following  result  is 
obtained:  4 

0  =  s  det(Is-A  +  BH)  +  k.^  N2-  (s)  h.-  (42) 

A  digital  program  is  used  to  solve  equation  (42)  for  the  closed 
loop  ACLS  system  roots.  The  value  of  k  is  then  determined  according  to 
the  deviation  of  ACLS  roots  from  the  basis  AMOAC  roots  and  system  response 
speed  of  the  ACLS  guidance  modes. 

This  value  of  k  is  then  incorporated  into  the  ACLS  simulation 
for  purely  demonstration  and  evaluation  purposes.  This  particular  guid¬ 
ance  loop  is  most  likely  not  the  best  available.  E"en  better  performance 
can  be  expected  with  a  more  complex  guidance  procedure. 

C.  Desired  Closed  Loop  System 

Physical  considerations  of  the  carrier  approach  problem  impose 
severe  constraints  on  the  location  of  the  closed  loop  denominator  and 
numerator  roots. 

1.  Phugoid  and  Short  Period  Roots 

In  the  final  phase  of  the  approach,  i.e.,  10  seconds  aft  of 
the  ramp,  the  airplane  should  be  capable  of  precise  response  to  altitude 
command  inputs  in  order  to  stay  in  phase  with' carrier  deck  motions.  This 
requirement  indicates  a  phugoid  bandwidth  of  a  magnitude  comparable  to 
the  center  frequency  of  the  carrier  deck  vertical  motion  spectrum  w 
For  adequate  altitude  control,  a  phugoid  damping  ratio  c5f  between  0.2  and 
0.4  is  required. 
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An  aircraft  carrier  moving  through  a  calm  sea  at  30  knots  sets 
up  vertical  gust  turbulence  with  a  low  frequency  spectrum.  The  carrier's 
pitching  and  heaving  motion  creates  additional  vertical  turbulence  of  a 
higher  frequency.  The  value  of  used  in  this  report  is  0.6  radians 
per  second,  the  average  of  the  values  found  in  references  (d)  and  (e) . 

The  vertical  gust  spectrum  data  were  taken  from  reference  (e)  which  in¬ 
dicates  vertical  turbulence  of  amplitude  1.75  feet  per  second  centered 
at  0.3  radians  per  second,  and  additional  turbulence  of  4.15  feet  per 
second  amplitude  centered  at  2.3  radians  per  second.  In  order  to  meet 
the  requirements  on  terminal  altitude  control  and  avoid  serious  vertical 
gust  problems,  the  phugoid  roots  should  lie  in  the  smaller  shaded  region 
of  figure  7.  Thus,  the  desired  phugoid  roots  are  characterized  by: 


0.2  S  t  5.  0  4 

0.47  S  oj  <•  0. 75 

P 

The  short  period  roots  should  have  a  bandwidth  smaller  than 


the  high  frequency  vertical  gusts  and  a  natural  frequency  great  enough 
to  provide  ample  separation  from  the  phugoid  roots,  i.e.,  wSp/Up  2.  2.5, 
For  satisfactory  pitch  control,  the  short  period  roots  should  have  a 
damping  ratio  of  between  0.55  and  0.85,  as  shown  in  the  larger  shaded 


region  in  figure  7  i.e.  desired  short  period  roots  are  char* cterized  by: 


0.5S  <  C  < 

0.85 
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2.  Nunerator  Zeros 

In  order  to  determine  the  significance  of  the  numerator  terms 
of  the  closed  loop  system,  it  is  assumed  that  pitch  rate  is  controlled 
primarily  by  the  elevator  deflection,  airspeed  primarily  by  thrust  and 
sink  rate  by  flaps  and  thrust  together.  Accordingly,  primary  interest 
is  given  to  those  transfer  functions. 

In  evaluating  the  zeros  of  the  sink  rate  and  airspeed  transfer  func¬ 
tions,  serious  attention  is  given  to  the  necessary  quick  glide  slope  re¬ 
sponse  and  the  non-minimum  phase  nature  of  the  open  loop  system  For 
these  reasons,  the  airspeed  and  sink  rate  transfer  functions  are  required 
to  have  all  left  half  plane  zeros.  In  addition,  the  magnitude  of  these 
zeros  was  limited  to  three  radians  per  second  in  order  not  to  lose  needed 
phase  lead.  The  pitch  rate  requirements  are  less  critical  and  therefore 
more  straightforward.  As  the  phugoid  and  short  period  roots  are  moved  to 
higher  frequencies,  the  pitch  rate  zeros  must  also  be  moved  to  higher 
frequencies.  This  procedure  insures  a  flatter  frequency  response  and 
makes  pitch  rate  response  less  sensitive  to  carrier  pitch  induced  dis¬ 
turbances.  This  phenomenon  is  illustrated  in  figure  8.  Thus,  the  mag¬ 
nitudes  of  the  two  pitch  rate  zeros  are  required  to  be  of  the  same  order 
as  the  natural  frequency  of  the  closed  loop  phugoid  and  shoit  period  roots 
respectively. 
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PART  III 
RESULTS 

A.  Application  to  the  F-8C 
1.  System  Synthesis 

The  effects  of  the  performance  index  parameters  upon  the  roots  of 
the  closed  loop  optimal  system  are  most  readily  shown  in  a  series  of  root 
locus  plots.  Figure  y  through  12  show  the  loci  of  the  roots  of  the  three 
inner  loops;  i.e.,  the  sink  rate,  pitch  rate,  and  cross  coupling  loops, 
as  functions  of  bi  and  b2  •  Root  loci  for  the  numerator  and  denominaror 
loops  appear  in  figures  13,  14,  and  IS.  The  loci  of  the  outer  loop  roots, 
or  the  final  closed  loop  airframe  roots  are  shown  in  figures  16,  17  and 
IS  Due  to  the  symmetry  of  the  root  square  locus  plots,  only  the  upper 
left-hand  quarter  of  the  s-planes  is  shown. 

In  the  carrier  approach,  the  primary  task  is  flight  path  control, 
with  pitch  attitude  control  as  an  important  secondary  consideration  In 
order  to  fulfill  the  primary  task,  considerable  attention  must  be  given  to 
shaping  the  phugoid  locus,  without  degrading  the  short  period  response 
The  rationale  behind  particular  loop  closures  is  to  position  the  phugoid 
roots  while  maintaining  satisfactory  short  period  damping  and  bandwidth, 
according  to  section  II  C. 

Inner  Loops 

*•  > 

The  three  inner  loops  have  bi  and  b2  as  common  parameters  and 
therefore  must  be  closed  simultaneously. 

Sink  Rate  Loop  (43a) 

A  -  A  *  A 

N22(s)  +  b;  N2i(s)  +  b2  N2 3 (s) 
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Cross  Coupling  Loop 

bi  kl  (s)  *  b2'  +  b,"  b2  ftfl(s)  (43b) 

Pitch  Kate  Loop 

A42  Cs)  +  bi  Ax*  j  (s)  +  b2  Ax, 3  (.s J  (45c) 

Substituting  values  from  Appendix  G,  the  inner  loops  become 
Sink  Rate  Loop 

4.8  x  10^a  *'-s  +  2.8 b>  <s2  +  0. 50s+0.51>  +  398b'  <s+0.23> 

1 

<S+3.82>  <s+3.40-/ +  5. j/b  <s+l).034>  ^s7  +  0.38+1.29* 

7 

Cross  Coupling,  Loop 

•  + 

1.51  x  1077  b  's(s+2.44)>  +  9.24  x  l(f9  b  <s(s^  0.12)>  + 

1  2 

>  > 

178.9  bi  l\->\s(s+0.04)> 

Pitch  Kate  Loop 

2.18  x  10~U<s(s+0.31)>  <s+0.89>  +  4.84  b  *s(s  +  0.52)> 

<S  +  0. 0G3>  +  1,27  X  10  3  b^  <s(s+3.18j'  •  s'U.'UO.V 
where  <-i,\s)>  denotes  F(s)  F(-s). 

Figure  9  shows  tne  loci  of  roots  for  the  sink  rate  loop  for  various 
*  * 

combinations  of  b^  and  b^.  Figure  10  is  an  enlargement  of  the  region  of 
interest  of  figure  9.  The  loci  of  the  cross  coupling  loop  roots  is  shown 
in  figure  11  and  figure  12  is  a  locus  of  the  pitch  rate  loop  rocts. 

The  sink  rate  and  cross  coupling  closed  loop  roots  form  the  open  loop 
poles  and  zeros,  respectively,  of  the  numerator  loop.  Pitch  rate  closed 
loop  roots  are  the  open  loop  zeros  for  the  denominator  loop.  The  cross 
coupling  loop  reflects  second  order  interactions  between  control  input 
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FIGURE  10.  COMPLEX  SINK  RATE  LOOP  LOCUS 
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pairs,  whereas  the  sink  rate  loop  reflects  a  primary  interaction  between 
sink  rate  and  the  control  inputs.  For  this  reason,  it  was  assumed  that 
the  cross  coupling  roots  would  affect  the  numerator  loop  much  less  than 
would  the  sink  rate  loop  roots. 

In  light  of  this  assumption,  the  closed  loop  sink  rate  roots  can  be 
considered  as  identical  to  the  closed  loop  numerator  roots.  Figure  9  shows 
the  entire  range  of  possible  roots  for  the  sink  rate  loop.  With  bi  and 
b2  equal  to  zero,  there  is  a  complex  root  pair  from  the  sink  rate  to 
thrust  transfer  function  called  phugoid  control  zeros,  and  a  negative  real 
root  at  s  =  2.86  rad/sec.  Since  the  phugoid  airframe  roots  exhibit  a 
small  natural  frequency,  a  complex  pair  of  zeros  will  be  required  to 
adequately  shape  the  phugoid  mode.  The  much  larger  natural  frequency 
of  the  short  period  mode  allows  for  satisfactory  shaping  with  a  single 
real  zero. 

Satisfactory  placing  of  the  phugoid  roots  must  include  a  large  increase 
in  speed  of  response,  increased  damping,  a  bandwidth  close  to  the  carrier 
deck  motion  center  frequency  and  adequate  separation  from  the  short  period 
roots.  These  physical  considerations  limit  the  complex  sink  rate  root  to 
the  shaded  region  of  figure  9.  Figure  10  is  an  enlarged  version  of  this 
desirable  region. 

Figure  12  shows  the  loci  of  the  pitch  rate  loop  roots.  For  all  values 
*  ■* 

of  bj  and  b2  there  are  three  negative  real  roots;  one  root  at  s  =  0, 

another  root  between  s  =  0.31  and  s  =  0.008  rad/sec  and  a  third  root 

>  > 

between  s  »  3.12  and  s  =  0.89  rad/sec.  Constraining  b]  and  b2  to  the 

same  values  as  in  figure  10  limits  the  two  negative  real  roots  to  the  indicated 

region  in  figure  12.  Using  any  values  of  bi  and  "t>2  in  this  limited  region 
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will  result  in  quite  similar  denominator  loop  root  loci. 

Three  sets  of  closed  loop  roots  were  chosen  to  illustrate  the  effect 

a  * 

of  large  ranges  of  bj  and  b?  on  subsequent  loop  closures.  Table  1  lists 
the  parameter  values  and  the  inner  loop  roots  for  these  three  cases. 


1 

CLOSED  LOOP  ROOTS 

Case 

a 

■ 

Sink  Rate 

Pitch  Rate 

Cross  Coupling 

I 

10' 13 

io'13 

-0.255  *j0.669 
-2.86 

0.0,  -0.305 
-0.885 

0.0,  -2.372 

II 

io'i2 

-0.342ij0.646 

-2.34 

0.0,  -0.293 
-0.844 

0.0,  -1.210 

III 

.  -12 
10 

0.362±j0.6SS 

-1.98 

0.0,  -0.287 
-0.860 

0.0,  -0-917 

TABLE  1  INNER  LOOP  CLOSURES 
Numerator  and  Denominator  Loops 

After  selecting  values  for  b1  and  b;>  it  is  now  possible  to  determine 
a  value  for  k2 . 

Numerator  Loop 

(Equation  (43a))  (Equation  (43b)  (44a j 

Denominator  Loop 

&  (s)  ♦  k  2  (Equation  (43c))  (44b) 

Substituting  values  from  Table  I  and  Appendix  G 
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Case  I 

Numerator  Loop 

8  .20  - 
4  808  X  Ilf  <s2+0. 510s+  5126>  <s+2.86.«-t-  1.602  X  10  k2  < s  (s+2. 372) ' 

Denominator  Loop 

<s2+0.96s+1.27>  <s2+0 . 012s+0.022>  +  2.23  X  10'nk2  s2  +  0. 305s>< s+0,885> 
Case  II 

Numerator  Loop 

8  _  1 9  ■* 

6-625  X  10  <s2  +0 ■ 683s  +  0-534>  *s  ♦  234>  +  6  2  X  10  k2  *-s(s+1.21)> 

Denominator  Loop 

<s2+0.96s  +  1.27>.^s2+0.012s+0.022>  +  2  670  X  10" ' J k2 ■ S  (s+0  844  )  -  <s+0-292- 
Case  III 

Numerator  Loop 

8<  2 

8.41  X  10"  s  +0-  724s  0.558^  <s*1.98'  ♦  1.09  X  10  1  ,Jk2  <s  (s+0.917)  ^ 
Denominator  Loop 

<s2+0.96s+1.27>  <s2+0.012s+0.Q22>  +  2.68  X  Ilf 1 ^s (s+0. 287) -  <s+0.860> 
Figures  13,  14  and  15  show  the  numerator  and  denominator  loci  for 
Cases  I,  II  and  III  respectively.  In  all  three  cases,  it  is  noted  that 
the  denominator  locus  experiences  much  greater  excursions  than  the  numer¬ 
ator  locus  for  equal  increments  in  k2. 

For  Cases  I  and  111,  the  open  and  closed  loop  numerator  roots  are 
nearly  identical  while  the  denominator  roots  experience  a  definite  ex¬ 
cursion,  i.e.,  the  phugoid  root  travels  approximately  15  percent  of  its 
total  locus.  This  evidence  validates  the  assumption  of  the  previous 
section;  that  is,  the  closed  loop  numerator  zeros  are  approximately  the 
sink  rate  zeros. 
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FIGURE  13.  NUMERATOR  AND  DENOMINATOR  LOCI  -  CASE  I 
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FIGURE  14.  NUMERATOR  AND  DENOMINATOR  LOCI  -  CASE  II 


Noting  the  limited  region  in  which  the  phugoid  locus  lies  and  relatively 
slow  movement  of  the  complex  numerator  root,  the  phugoid  locus  in  the 
outer  loop  closure  can  be  expected  to  be  largely  independent  of  the  kji 
closure.  The  short  period  locus  displays  markedly  different  behavior  than 
the  phugoid  locus.  Evidently,  values  of  k2  can  have  a  considerable  effect 
on  the  short  period  !icus  in  the  outer  loop  closure.  In  order  to  demon¬ 
strate  this  conjecture,  Cases  1  and  III  are  closed  around  small  values  of 
k^  and  Case  II  is  closed  about  a  larger  value  of  k',  as  shown  in  Table  II. 


CLOSED  LOOP  ROOTS 

B 

_ 

Numerator 

Denominator 

i 

o 

o 

H 

-0.256  *  j0.668 
-2.86 

-0.053  *  jO. 134 
-0.573  *  j 1.019 

a 

xou 

-0.419  *  j0.574 
-2.47 

-0.114  1  j0.038 
-1.105  *  j 0.857 

hi 

►-* 

o 

e— 

O 

-0.372  *  jO. 644 
-2.00 

-0.053  *  jO. 133 
-0.587  *  j 1.016 

TABLE  II.  NUMERATOR  AND  DENOMINATOR  LOOP  CLOSURES 


The  closed  loop  roots  can  now  be  specified  with  a  single  loop  closure. 
OUTER  LOOP 

(Equation  44b)  +  k{  (Equation  44a)  (45) 

Substituting  values  from  Table  II 
Case  I 


<s2+0 . 106s* . Q205>  <s2+l. 146s+l. 367>  ♦  4.81  x  10"Bkf<s+2.86>  <s2*0.512s*0.512:> 


Case  II 

<s20.228s+.0146>  <s2+2.21s+1.95S>  +  6.63  x  10_8ki  <s+2.47>  ^s2  +  0.838s*0. 515> 
Case  111 

<S2+U.  106s+.02Q3>  <-s2  +1 . 164s+l .  377>  +  8.42  x  10'8  ki  <s+2.0>  <s2+0- 744s+0. 553' 

Figures  16,  17  and  18  show  the  loci  of  the  closed  loop  roots  for  Cases  I, 

II  and  III  respectively,  The  closed  loop  short  period  and  phugoid  roots 

a 

are  shown  for  various  values  of  ki .  For  convenience  of  comparison.  Table 

•  *  a 

III  lists  these  roots  and  their  corresponding  values  of  bi ,  b2,  ki  and 

h- 

The  phugoid  locus  of  the  outer  loop  closure  is  most  greatly  affected 
by  the  complex  closed  loop  numerator  zero.  This  zero  is  essentially  the 
same  as  the  closed  loop  sink  rate  zero.  Thus,  the  closed  loop  phugoid 
locus  can  be  shaped  effectively  by  the  sink  rate  loop  closure.  The  short 
period  roots  of  the  outer  loop  move  considerably  more  slowly  along  their 
locus  than  do  the  ph'^oid  roots.  This  places  greater  emphasis  on  the  ini¬ 
tial  part  of  the  short  period  locus  than  upon  latter  parts.  Therefore,  the 
short  period  locus  is  greatly  affected  by  the  closed  loop  denominator  roots. 

The  design  procedure  used  in  this  report  allows  a  great  deal  of  independ¬ 
ence  in  placing  the  closed  loop  short  period  and  phugoid  roots.  The  closed 
loop  phugoid  locus  is  shaped  primarily  by  the  sink  rate  loop  closure,  i.e. 
the  values  of  bj  and  b2,  while  the  closed  loop  short  period  locus  is  in¬ 
fluenced  more  by  the  denomin? tor  loop  closure,  i.e.  the  value  of  k2. 


FTf.URE  18.  FINAL  LOOP  CLOSURE 
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2 .  Design  Evaluation  , 

In  this  section,  .the  preliminary  design  of  the  previous  section  is 

.  i 

evaluated  ahd  modified  according  to  the  procedure  in  section  II  in  order  to 
achieve  the  design  of  a  fully  automatic  approach  controller  for  the  longi¬ 
tudinal  axis  of  an  F- SC  aircraft  In  the  following  discussion  various 
candidate  design  configurations  will  be1  introduced.  In  order  to  minimi ze 
confusion,  each  of  the  preliminary  designs  will  be  designated  by  a  system 
number  This  number  is  a  chronological  ordering  of  the  candidate  designs. 
The  baseline  system  is  number  zero  with  six  other  candidates  systems  com¬ 
prising  the  entire  list  1 

1 

The  results  of  the  previous  section  indicates  a  baseline  system 
suitable  for  initializing  theiAMOAC  design  procedure.  The  particular  system 
is  Item  7  of  table  III.  For  convenience,  table  IV  below,  contains  a  list 
of  the  seven- preliminary  designs  along  with  their  respective  performance 
index  parameters. 


System 

No. 

am 

n 

* 

ki 

- 

1  '  ■ 
k2 

0 

io‘n 

lu’13 

— 

4 

X 

107 

i.jio 

j 

1 

10"  1 1 

10"  1? 

4 

X 

IQ7 

• 

idi0 

2 

id"  n 

2  x  K>';2 

4 

X 

107 

l 

to'0  , 

3 

10" ir 

10*  li 

4 

X 

10; 

10' 0 

4 

10' 13 

2  x  10_u 

4 

X 

107 

io10  j 

5 

10' 1  * 

10" 

4 

X 

107 

10'iu 

6 

10'  13 

2  x  10‘!0 

4 

X 

107 

10 !0 

TABLE  IV  -  PRELIMINARY  SYSTEMS 
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For  eaUi  system  listed  Ln  table  IV,  a  feedback  gain  mutrix  and  a  closed 
loop  transfer  matrix  lias  been  calculated  as  described  in  section  I  C. 

These  data  are  listed  in  Appendix  H. 

Inspection  of  the  transfer  matrix  of  system  0  reveals  that  the  sink 

rate  to  thrust  transfer  function  has  right  half  plane  zeros.  However, 

the  other  requirements  were  fulfilled.  These  non-minimum  phase  zeros  are 

a  complex  pair,  close  to  the  imaginary  axis,  close  in  magnitude  to  the 

left  half  plane  sink  rate  to  flap  zeros.  For  this  reason,  it  was  felt  that 

a  greater  use  of  flaps  as  opposed  to  thrust  would  result  in  left  half  plane 

sink  rate  zeros  To  accomplish  this,  the  performance  index  weighting  of 

flap  deflection  was  decreased  while  the  other  parameters  remained  unchanged 

** 

The  parameter  b2  the  reciprocal  of  d.,  was  varied  from  the  system  0 
value  of  10~13  through  four  orders  of  magnitude  to  a  value  of  10_i.  The 
root  square  locus  was  used  to  monitor  movement  of  the  closed  loop  system 
roots.  This  root  locus  plot  is  shown  in  figure  19. 

A 

Six  different  values  of  were  chosen  to  define  the  six  candidate 
designs  designated  system  1  tnrough  system  6.  Inspection  of  the  transfer 
matrices  of  these  closed  ioop  systems  eliminated  consideration  of  candidate 
systems  1,  2,  3  and  4  since  each  of  these  systems  contained  right  half 
plane  zeros  in  their  sink  rate  to  thrust  transfer  functions  System  S  and 
system  0,  however,  satisfied  the  criteria  specified  for  satisfactory  nu¬ 
merator  roots. 

The  next  design  step  is  the  exercise  of  a  carrier  approach,  simulating 
the  simplified  aircraft  model,  with  AMOAC  augmentation  by  each  of  the  two 
candidate  systems  respectively,  as  described  in  section  II  B  -2.  Figures 
20  and  2i  contain  the  time  histories  of  the  carrier  approach  employing 
system  S  and  system  6  respectively.  Inspection  of  these  time  histories 
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FieURE(20r  SYSTEM  5  CARRIER  APPROACH 


shows  a  generally  acceptable  response.  The  glide  slope  angle  shows  no 
appreciable  degradation  up  until  one  hundred  feet  to  touchdown  and  the 
variations  from  trim  velocities  remain  tightly  bounded,  as  compared  to 
the  magnitude  of  the  disturbances.  The  height  over  ramp  and  touchdown 
dispersion  criteria  were  satisfied,  as  shown  in  table  V. 


Pitch  rate  response  and  elevator  deflection  were  minimal  The  DLC 
flap  deflection  remained  within  operational  bounds;  however,  the  thrust 
limitations  were  exceeded.  Since  the  periods  of  thrust  limit  exceedance 
were  of  short  duration,  the  response  was  classed  as  acceptable  and  the 
next  step  in  the  design  procedure  was  initiated. 

The  airplane  modei  was  expanded  tO' include  lags  from  the  control  sur 
face  actuators  and  engine,  and  the  hard  limiters  on  the  controller  output 
The  specific  values  of  the  lags  and  limits  are  included  in  appendix  D 
Figure  22  and  23  show  time  histories  for  systems  5  and  6  respectively. 
These  time  histories  are  quite  similar  to  the  previous  set.  The  criteria 
on  AHp  and 

ity  variations  are  again  well  bounded  and  in  general  smoother  than  the 
previous  set.  Thrust  output  is  only  limited  over  a  small  segment  of  the 


are  still  satisfied,  as  shown  in  table  VI.  The  fiin  velcc 
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response  and  exhibits  linear  control  over  the  majority  of  the  final  one 
thousand  feet.  The  DLC  flap  controller  operates  entirely  in  its  linear 
region.  With  these  simulations  satisfactorily  completed,  the  two  candi- 


Equation  (41)  defines  the  simple  ACLS  guidance  equation  used  in  the 
ACLS  simulation.  In  order  to  find  a  proper  value  of  guidance  gain  k,  a 
root  locus  of  the  guidance  loop  was  constructed  according  to  equation  (42). 
Figures  24  and  25  show  the  root  loci  for  systems  5  and  6  respectively. 

In  order  to  be  consistent  with  the  ten  second  recovery  requirement  on 
glide  slope  control,  closed  loop  system  roots  should  have  negative  real 
parts  with  magnitudes  greater  than  0.2  radians  per  second;  however, 
magnitudes  in  excess  of  1.0  radian  per  second  seem  unnecessary  Inspection 
of  figures  24  and  25  show  that  only  guidance  gain  values  between  0,2  and 
0.6  need  be  investigated.  Figures  26,  27  and  28  show  20  feet  low  approach^ 
of  system  5  with  guidance  gains  of  0.2,  0.25  and  0.3  respectively.  Figures 
29,  30,  21  and  32  show  20  feet  low  approaches  of  system  6  for  guidance 
gains  of  0.2,  0.25,  OcoO.and  0.35  respectively. 

Inspection  of  figures  26  through  32  show  the  expected  quickening  of 
glide  slope  recovery  with  increasing  guidance  gain  k.  However,  as  k  is 
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o.o  —  tiior  icwr  8Mtcr  toco) 
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increased,  non-linear  effects  build  up  and  prolong  the  recovery  of  airspeed 
and  sink  rate.  In  particular,  the  large  assymetry  in  thrust  limits  com¬ 
bines  with  the  large  engine  lag  effect  airspeed  response  adversely.  As 
k  is  increased,  the  system  uses  the  large  positive  thrust  limit  to  null 
positional  errors  quickly.  The  resulting  speed-up  takes  considerable 
time  to  null  out  when  thrust  is  reduced  to  its  small  negative  limit- 
Since  the  glide  slope  response  time  and  the  airspeed  response  tim  trends 
oppose  each  other,  there  must  e:;ist  a  value  of  k  for  each  system  with 
fastest  overall  response.  Figure  33  shows  these  trends  and  the  resulting 
optimum  response  times  Inspection  of  figure  33  shows  that  system  6  re¬ 
sponds  1.6  seconds  faster  than  system  5.  Though  neither  system  can  match 
the  ten  second  criterion,  figures  34  and  3S  show  that  system  6  comes  with¬ 
in  0-2  seconds  of  compliance,  and  engineering  judgment  dictates  acceptance 
of  this  performance.  Thus,  sy sr  5  was  dropped  from  consideration  for 
failing  to  provide  fast  glide  slope  response,  and  system  6  was  retained. 

The  ACLS  employing  system  6  and  the  optimum  guidance  gam,  k  *  0.24,  was 
then  simulated  in  a  fully  automatic  carrier  approach.  Figure  36  shows  a 
carxier  approach  from  a  20  foot  initial  glide  slope  error,  through  the 
burble  to  the  TOP.  For  this  run,  all  of  the  landing  requirements  were 
satisfied.  The  height  over  tne  ramp  error,  AIL,  was  3.8  feet  and  the 
touchdown  height  error,  was  2.34  feet.  The  carrier  motion  height 

dispersion  was  calculated  according  to  equation  (39)  to  be  0.2  feet, 
a  small  number  relative  to  the  design  criterion. 

Ihus,  tiie  control  law  represented  by  the  system  b  feedback  matrix  has 
been  shown  to  be  an  effective  carrier  approach  controller,  and  the  AMOAC 
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e.O  _  CLIDE  M.BPE  MULE  10ECJ 
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PART  IV 
CONCLUSIONS 

The  feasibility  of  the  AMOAC  design  procedure  has  been  demonstrated 
by  its  application  to  the  control  of  the  longitudinal  axis  of  the  F-8C 
fighter  aircraft. 

The  blending  of  optimal  control  theory  and  modern  state  variable 
analysis  techniques  have  allowed  the  synthesis  of  a  complex  twelve  para¬ 
meter  feedback  control  system,  with  the  ease  and  simplicity  of  a  four 
parameter  design  technique.  The  iterative  nature  of  the  design  logic  was 
shown  to  converge  rapidly  in  the  F-8C  application.  After  the  original  use 
of  the  root  square  locus  this  design  tool  was  used  only  once  in  updating 
the  design  parameters.  In  spite  of  the  imp’ icit  assumption  of  system  line¬ 
arity,  the  AMUAC  design  procedure  was  successful  in  treating  the  non- 
linearities  of  the  F-8C  airframe  and  propulsion.  The  resulting  glide  dope 
control  is* very  precise.  Variations  of  aircraft  position  from  the  nominal 
glide  slope  is  held  to  a  tight  tolerance  in  the  presence  of  severe  carrier 
burble.  The  AMOAC  control  system  has  demonstrated  its  ability  to  correct 
a  20  foot  vertical  glide  slope  error  in  under  eleven  seconds.  In  addition, 
this  system  responds  well  to  carrier  deck  motions.  The  use  of  an  AMOAC 
system  in  actual  carrier  approaches  should  yield  a  substantial  improvement 
to  present  automatic  carrier  landing  capability  The  results  o I  Appendix  I 
show  tlvat  similar  closed  loop  dynamics  can  be  expected  for  other  carrier 
based,  high  performance  aircraft.  In  particular,  it  is  fe'lt  that  the 
AMOAC  design  procedure  has  potential  for  rapid  and  straight  .iotward  ap¬ 
plication  to  the  A-5A  and  the  F-4  aircraft. 
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The  results  indicate  that  the  following  areas  should  be  investigated 
in  greater  depth-  These  areas  are:  determining  the  effects  of  carrier- 
deck  motion  prediction,  determining  the  effects  of  deleting  the  less  im¬ 
portant  feedback  gains  and  the  applying  an  AMOAC  type  design  procedure  to 
the  1 iteral-directional  mode  of  a  powered  approach, 
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AERODYNAMIC  PARAMETERS 


A  Definition  of  Dimensional  Derivative! 


'u  mVo  + 


S  q 
wn 


*w  =  — - - (C,  -  cn  ) 


mVo 


L  D„ 


X  =  V0  X 
a  w 


S  q 

v  -  w 
6li  m 


06,: 


X  t  =  --£*a -  Cn 


6T 


£!£l_iiw  +  0 
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mV 
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(CL  +  V 
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mV0. 
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am 
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sin  (iw  ) 
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B.  Aerodynamics  of  Power  Approach  (F-8C) 

Physical  Constants 

W  =  20,000  lbs 

Iy  =  90,000  slug- ft2 

S  =  375  ft2 

w 

c  =  11.78  ft 

iw  -  7.0  deg 

K  =  0.85  dep 

l  =  -0.424  ft 

Initial  Conditions 

V0  =  134  knots  =  226  ft/soc 

Y0  =  -4  deg 

«o  =  =  13.6  deg 

“fus  a  \  \  •  6-6  de« 

9o  =  ao  +  Yo  =  9-6  deg 

5E0  =  -7  deg 

5F0  =  -15  deg 

To  =  2494  lbs  e  80%  rpm 

C.  Dimensional  Derivatives  for  F-8C 

For  q  =  60.7  lbs/ft.,  and  c.g.  =  27  percent  of  c. 

Xu  -0-0506 

X  =  -0.0072 

w 

X  =  1.63 

a 

X6E  -  -0.672 

XfiF  =  -6.72 
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X6T  =  0.001S95 

Zu  s  Q.264 

Z  =  -0.554 

w 

ZXL.  =  -19.95 

o  c 

Z6F  *  '18-9 


Zx...  =  -0.0002194 

6  1 

M  B  1.074 

a 

M  *  -0.000189 

w 

M  =  -0.043 

a 

-  -0.347 

M...  =  -2.20 

o  c 

M6F  =  ° 

M6T  +  -4.71  x  !0'6 


0 .  Approach  Parameters  for  Various  High  Performance  Aircraft 
Physical  Constanta 


Airplane 
I  (.lbs) 
ly  (slug-ft  ) 

sw  (ft2) 

(ft) 
iw  (deg) 

4  (deg) 
lz  (ft) 
c.g  (*MAC) 

V*  (ft/sec) 


1*  -  4 11 

A-5A 

F-111B 

34,000 

38,500 

57,800 

111,273 

252,512 

312,000 

530 

700 

550 

16.04 

15.2 

8.8 

1.0 

0 

1.0 

4.25 

1.24 

1.20 

-0.364 

0.08 

0.71 

30.0 

30.6 

35.5 

226 

223 

200 
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Dimensional  Derivatives 


A i rplane 

F-4B 

A-SA 

X 

u 

-0.094 

-0-0544 

-0 

X 

a 

4.94 

12.57 

16 

X<ST 

9.43  X  10~4 

8.36  x  10"4 

5 

X6F 

N,A. 

-5,16 

-4 

z 

u 

-0.305 

-0.274 

-0 

L 

a 

-83.6 

-146.7 

-99 

Vi 

-8.67  x  10"5 

-1,82  x  10-5 

-1. 

“6F 

N.A. 

-26.75 

-32. 

1  u 

1.8  x  10-4 

-8.0  x  10"5 

-2- 

-1.298 

-0,88 

-0. 

-3.49  x  10"4 

-4.62  x  10"4 

-4. 

1 

a 

-0.0789 

-0.103 

-o. 

1 

q 

-0.433 

-0.288 

-0. 

l6T 

-3.27  x  10-6 

*3.7  x  10"8 

2. 

F-nib 

.0456 

.45 

.57  x  10 
.93 

,309 

.8 

.89  x  10' 

25 

05  x  10' 
08 

28  x  10" 

0856 

436 

26  x  10  ‘ 


OTE:  N,A.  means  Not  Applicable 
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APPENDIX  B 


A.  Euler  Equations 

For  a  given  set  of  differential  constraints  of  the  form 

x  =  f  (x,u)  x(0)  =  c  (B-l) 

and  a  convex  functional,  J  (jc,u)  ,  of  the  form 

J(x,u)  =  Jo  L  (x,u)  dt  (B-2) 

it  is  desired  to  find  vector  functions  u*(t)  and  £*(t)  that  simultaneously 
satisfy  equation  (B-l)  and  minimize  J(x,u).  For  this  development,  the 
constraint  equation  (B  -1)  becomes  the  open  loop,  longitudinal  airframe 
equation  evaluated  along  the  nominal  approach  glide  path.  The  vector  x^ 
is  expressed  in  the  airplane’s  stability  axis  coordinates,  see  figure  (1). 
The  functional  J(x,u),  is  the  performance  index  defined  by  equation  (4). 


x  *  Ax  +  B6 


/■  oa 

J  =  1/2  J  [X  fK  X  +  61  DS  ]  dt 


x(0)  sc  (B-3a) 
(B- 3b) 


A  new  set  of  functions,  £  (t)  are  introduced  to  adjoin  the  constraint 


equation  to  the  performance  index. 

J'  =  1/ 2  [x1  Kx  +  £lD£  -  2£F  (x*  -Ax  -  B£)  ]  dt  (B-4) 

It  should  be  noted  that  when  x^  and  £  satisfy  (B-3a),  J  =  J  regardless 
of  £.  The  calculus  of  variations  states  that  a  necessary  condition  for  J 
to  take  on  a  minimum  value  is  that  the  first  order  difference  of  J  (x*'» 

&*)  J(x,£)  be  zero 

A  J  *  J  (x*,  6*)  -  J(x,  6) 

°  1/2  I  o  [x*TKx*  -  x1  Kx  ♦  £*TD6*  -  _6T D6_  + 


*2E  (x*  -  x  -  Ax*  ■»  Ax  -  B6*  B£)  ]  dt 


(B-S) 
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Using  Taylor's  thee  em 

T  T  T 

x*  Kx*  -  x  Kx  =  (2Kx*1  (x*  -  jc)  +  higher  order  terms 


=  (2Kx*)r  Ax  ♦  0  C  I A x  |  ^  ■) 


Si  ii  iarly 

A^OS*  -  iT  *  (2D£*1)A6^  +  0(lA6|?) 


Substituting  into  equation  (B-31 

,  -*  T  T 

AJ  =  \  [  (Kx*  +  A1  P)  AX  +  (D6.*  +  B*  P)  A6  -  pTAx  ]  dt 

* j o  CO  (  I  Axj  ? )  +  Of  I A6 |2)]dt 


Using  integration  by  parts 

r  •'o  ^  Too 

I  £*  Ax  dt  =  -  I  Ax  dt 


T 

♦  £  A  x_ 


GO 


o 


Since  there  is  no  perturbation  in  the  initial  value  of  £(t),  i.e  ,  Ax(0,i 

=  (0)  equation  (B-S)  becomes 

C  o<  .  T  ,  t 

AJ  r  [  (Kx*  +  A!£  +  £)  Ax  +  (D6*  +  BrPJ  A6j  dt 

£J  («]  Ax.  (”)  +  J  [0  ( |  Ax_| 2)  ♦  0(|A6|*)]  dt  (l'-bj 

The  relation  between  ax  and  A6_  is 


Ax_  =  aA£  +  BA£  Ax(0)  =  0 

SinceA£  is  unconstrained  but  small  in  magnitude  the  resulting  Ax^ 
spans  the  entire  x  region  about  x/ .  Thus,  the  requirement  that  *he  fjrsi 
order  difference  vanish  yields  the  following  necessary  cond’ t  w.  on  y_" 
and  6*. 


£  =  -  AT  £  -  K  X* 

x(“)  =  0 

(B-.7a) 

0  =  D6*  +  BT£ 

.  (B* 7b) 

x  =  Ax  *  +  B6* 

x(0)  =>  c 

(B-7c) 
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Equations  (B-7a,  b  and  c)  are  the  Euler  equations  for  the  optimal  system. 
B.  Riccati  Equation  and  Feedback  Gains 

The  solution  of  equations  (B-7a,  b  and  c)  has  the  form 


(t) 


x*(t) 


*11  Ct) 


*21  (t) 


*12(0 


*22  Ct) 


2.(0) 

x(0) 

X  (“)  =0 


Solving  this  equation  for  p  in  terms  of  x  yields 

£(t)  =  -  [*ll(-t)  +  r  *2l(-t)]'1  [  *  1 2  ( "  t )  ♦  ( *22  C*t)J  X*(t) 

where  T  =  lim  *21  (-T)  *22  (TJ 

T-*® 


Thus 

£(t)  =  p(t)  x*(t)  (B-8) 

where  P(t)  is  the  n  x  n  Ricatti  matrix.  Combining  equation  (B-7c)  and 
equation  (B-7b) 

£  =  -  Ar  £  -  Kx* 


Ax*  -  B  D_1Br£ 


(B-9o) 

X"  =  Ax’*  -  B  U  p  (B-9b) 

Combining  equations  (B-8)  and  (B-9)  gives 

-Px*  =  (PA  +  AT  P  +  K  -  PBlf1  B!P)  x*  (B-10) 

Since  equation  (B-10)  is  independent  of  the  initial  choice  of  x,  the 
Riccati  matrix  must  satisfy  the  following  equation. 

-P  =  PA  +  AT  P  +  K  -  PBD-1  Brp  (B-ll) 

Equation  (B-7b)  shows  the  form  of  the  optimal  feedback  control  law 

(B-12) 


£*(t)  =  -O"1  Bf  £  =  -D'1  BTP  (t)x  (t) 
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For  this  report,  the  time  interval  of  interest  is  0  t  'S  ns  such, 
the  optimal  control  6_ * ( t }  corresponding  to  a  disturbance  x(t)  and  the 
optimal  c 'Pt  1  ■  ■  iS?  fl  corresponding  to  x  (t  *  <0  are  related  by  £2*(t) 

■»  a  *  ft  •  ')  (B-13) 

This  states  the  feedback  control  law  is  time  invariant  or  that  the 
Riccari  matrix  is  a  matrix  of  constants  Thus  P ft  1  =  0  and  the  Riccati 
equation  reduces  to 

0  -  PA  +  AT  P  *  K  -  PBD*  BPP  (B-I4- 

and 

6*(tl  =  -D~  BfP  X  35> 
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APPENDIX  C 

ROOT  SQUARE  LOCUS  EQUATION 


In  Appendix  B,  the  differential  equations  of  the  optimal  system, 
equations  (C-7a,  b  and  c)  were  derived,  and  may  be  written  as: 


ABO 

T 

0  D  B 

T 

-K  0  -A 


(C-l) 


This  equation  completely  determines  the  closed  loop  roots  of  the 
control  system  and  its  adjoint.  In  order  to  demonstrate  the  nature  of 
the  closed  loop  roots,  it  is  necessary  to  solve  the  characteristic  poly¬ 
nomial  of  equation  (C-l).  Laplace  transforming  equation  (C-l)  yields: 

Is  -  A  -B  0 


0 
-  K 


-B 

-D 

0 


•B 


-IS 


.T 


-] 

x(s) 

£00 

-£(s) 

.  *0 


(C-2) 


whose  characteristic  polynomial  is: 

Is  -A  -B 

0  -D 

-K  0 


det 


=  0 


(C-3) 


This  (2n  ♦  ■)  dimensional  matrix  is  divided  into  nine  parts.  The  first 
n  rows  correspond  to  the  airframe;  the  next  m  rows  to  the  control  vector 
A(s),  and  the  last  n  rows  to  the  adjoint  variables  £  (s).  Adding  a 
linear  combination  of  the  first  n  rows  to  the  last  n  ro. '  does  not  alter 
the  determinant.  Thus  equation  (C-3)  can  be  simplified  to: 
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dot  (D  ♦  M(-s)TKM(s)j  *  0 

Equation  (C-61  is  the  basic  equation  of  the  root  square  locus 


(C-6) 


'  —  “1 
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APPENDIX  D 

;  \  , 

Carrier  Approach  Simulation  Equations 

For  definitions  of  positicfa  coordinates  and  axis  systems,  refer  to 

figure  (6). 

•  I  *  1 

A.  Positional,  Coordinates 

!  i  ; 

The  nominal  glide  slope  i,s  a  four  degreo  glide  path  starting  6000.0 

‘  ■  1  •  ’  , 

feet  behind  and  419.4  feet  above  the  touchdown  point  of  a  stationary  aii- 

I 

craft  carrier.  The  aircraft's  closing  velocity  is  180.0  fer-r  per  second 

1 

Upper  case  letters  denote  nominal  ghde  slope  quantities  and  I'ower 

»  t  , 

;  i 

case  letters  are  aircraft  centered  variables. 

» 

X  -  6000. 0  -  180. Ot  .  (D-^a) 

Z  =  419.4  -  12. 58t '  (D- lb) 

The  aircraft's  deviation  from  the  glide  slope  ij  i 


Ax 

u  dt 

(D*-2a ) 

ft 

1 

.  '  i  i  ‘ 

!  ■  i 

Ah  x  J  o 

h  dt  i 

,  (D-2b) 

Relative  to  the  carrier,  the  aircraft  is  positioned  as  follows: 

*  *  X  ♦  Ax  •  ,  «  fD-*3a) 

-  h  .*  Z  ♦  Ah  C0-3b) 

R  *  (x;+h2)2.  (D-  3c  V 

I 

T  ^  tan"1  (h/x)  ,  (D-3d) > 

B.  Aircraft  State  Variables 

The  state  variables1 of  the  aircraft  are  the  same  as  in  Section  U  A. 

*  I 

At  each  instant  of  time  the  starred  quantities  represent  the  aircraft 
variables  relative  to  the  aircraft  carrier  burble. 


1 
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where* 


u*  =  u  ♦  Ug(R) 
h  ■■  h  *  Wg(R) 


-R/1154 


23.0eK/1154-  17.5e«/75 


(1.8  -  0.0Q9R)  sin  (26.0-R/41 . 7.) 


(D*  4a) 
(U-4b) 

R  ^2000  i't 
R  *■  2000  ft 


-7.0  eR/ll54cos  (R/255)  R  .i.2306  ft 
-7.0  eR/llS4cos  (R/255)  R  r  2306  ft 
♦(4. 15-0.00185R)  sin  (26.0  -  R/41.7) 


The  aircraft's  equations  of  motion  are  then 

x  =  A  x’ »  ♦  B  ^  .  (D-5  ) 

C.  Control  Law 

When  the  aircraft  is  not  engaged  in  a  closed  loop  ACLS,  the  control 
law  is  identical  to  the  one  in  equation  (8).  Thus: 

i  r  -ii  x*  rn-6) 

is  used  in  the  simulation. 

When  the  ACLS  mode  is  used,  a  quantity  proportional  to  the  vertical  glide 
slope  error  is  added  to  the  sink  rate  feedback  term. 

s  h*  +  k (h- 0 . 069U )  10-7) 

Thus  the  ACLS  control  law  is: 

6  -  -H  xfb  (D-8.) 

D.  Airframe  Model 

The  linearized  airframe  equation  (U-5)  is  employed  only  in  *■  initial 
simulation  exercise.  Further  exercised  use  a  more  complicated,  nonlinear 
model.  This  model  includes  first  order  lags  for  elevator  and  flap  actu¬ 
ators  and  engine  thrust  build  up,  and  control  authority  limiters- 
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In  the  following  equations,  upper  case  subscripts  denote  controller 
command  and  lower  case  subscripts  denote  controller  outputs,  primed  quan¬ 
tities  have  been  subjected  to  authority  limitations 


Actuators 


5  '■  -L  ♦  l*_j  (0-10; 


Limiters 


I 

L  0 


0 


0 


(Dll) 


where  figure  (D-ll  illustrates  the  functional  relationship  Ae(  ),  Atf  ) 


K  I- 
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APPENDIX  E 

Carrier  Motion  Dispersion 

Appendix  3  described  the  system  equations  of  the  AMOAC  controller 

used  in  the  ACLS  mode.  In  this  appendix  a  method  is  derived  for  comput- 

•x. 

ing  the  carrier  motion  induced  glide  slope  error.  If  x  is  a  vector  of 
the  aircraft's  state  variables  plus  position  variables 


x_  =  column  (x,h,u,h,0,G) 

and  x  is  a  vector  of  the  aircraft  carrier's  motion 
-cm 


column  (x  ,h  ,0,0,0  01 

v  cm  cm*  *  an  * 


(E- 1) 


(E-2) 


and  is  a  unit  vector  with  ceros  in  .ill  but.  the  i-th  position,  i.e. 


ej  =  column  (0,1,0,0,0,01 

The;;  the  desired  vertical  error  term  is 

h  .  =  h-h  -  e1  (x  -x  ) 
d  cm  —  —cm 


(E-3j 


(E*  4) 


A  good  measure  of  vertical  dispersion  is  the  one  sigma  variable  distance 
’  f”  1  f 

-  j  ShJ  w  d»  (U-S) 


''  lid 
where 


r 

r 

1  dl 

lim  1 

/-OO 

T-  ;>T 

no  v 

T  -  .  ox  ; 

3hdvu'-'  J  u‘  |  ~  |  hdft41)hd(t-)dt  e  j 

_  -T 

is  the  power  spectral  density  of  the  glide  *  lope  error. 


(E-6) 


The  airframe  equations  are  as  follows: 

>.  v  > 

x  -  Ax  *  B  o 


'V  ■'  t  'C 

_6  =  -H(x-x) 


(Li-6,1 

(E-7) 


Combining  equations  (E-6)  and  (E-7)  and  Fourier  transforming  yields 

.  *v.  -V  _  Vt-V, 

X  *  (Ij  -A  ♦  BH)  BHX 


(E-8) 
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Laplace  transforming  equation  (E-4)  and  substituting  equation  (E-8)  yields 
Hd(to)  =  T  [(la>-A  *  BH) " 1  BH-ljx  (E-9) 

Upon  rearranging 

Hj  (u)  =  ejA(w)  k  foil  (E- 10) 

where 

A  (u»)  =  [f,wI-A  +  BH)  ” 1  (,  tuI-A)  ] 

Expanding  equation  (E-10)  yields 

Hd(w)  -  A ^  i  (to )  xcnj(w)  +  Hcm  (w)  «  A?5(u))0cm((j>) 

where 


A  (w)  =  e^  A(u>)  ej 


ij 

Thus  the  resulting  power  spectral  density  is 

2 


shd‘“> 


(A21  (u) 
A,>  • « to) 


W*>4 


(E-12 ) 


A2t>fu)) 


Sgcm(oi.)  ♦  Cross  Spectral  density 
terms 


Since  the  complete  set  of  cross  spectral  density  functions  is  not 
available,  their  contribution  was  dropped  from  equation  (E-12),  thus 
yielding: 


hd  " 


/I 

% 

CO 

- 

? 

( 

A2j (w) 

S  («'  * 
xcm  - 

A 22  (<*>) 

VJ 

oc* 

L 

SHcm  * 


A  2  i  (<*■') 


S0cm^ 


du>  2 


(E- 133 


All  the  matricies  and  functions  used  in  obtaining  equation  (E-13'» 


are  defined  in  section  V  of  this  report. 


pij  hi  r  ,i  -, 
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Appendix  F 

Guidance  Loop  Equations 

When  the  AMOAC  control  system  is  engaged  in  an  ACIS  mode,  the  system 
behavior  is  altered  m  accordance  with  the  guidance  loop.  An  equation  for 
the  closed  loop  ACLS  roots  is  derived  in  this  section. 

The  optimal  system  is  characterized  by 

x  --  A  x  ■-  B  5  (F  i.al 

--  -H  x.  (Fib} 

where 

x  -  column  (u,h,0,0)  and  6_  -  column  <„} 

The  guidance  command  adds  vertical  position  error  to  the  sink  rate  feedba; 1 . 


h.  *  h  ♦  k  h 
fb 


fp-2) 


Laplace  transforming  equation  (F-2)  and  combining  with  equation  (Fib)  yields 
u  =  H Cl  ♦  e  eT  £)  x  (1-3) 

where 

e  *  column  (0,1 .0,0) . 

And  the  closed  loop  ACLS  equation  becomes 

£ls-A  ♦  BH(I  ♦  e?ej  x  -  0  (F  41 

Upon  rearranging  equation  (F  4) 

1^  £  (Is-A  ♦  BH)s  ♦  kBH  e?e?j  x  =  0  (F-S.: 

Equation  (F-5)  is  obtained  which  results  in  the  characteristic  equation. 


det 


(Is-A  ♦  BH)s  ♦  kBH  e  er  = 

Lv  -2—?  J 


0 


or 


det  (Is-A  ♦  BH)  det  (Is  ♦  k(Js-A  +  BH)~  BH  *>,e?)  =  0 
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CF-7 : 


fe 


NADC-SD-7153 


However,  the  first  fa>  tor  in  equation  (F-7)  is  the  open  loop  charae 
terirtic  equation  and  equals  zero  only  when  <  equals  zero,  thus 


Jet 


[is  -  k  (Is-A  ♦  BH) *  * B  J  H  j  r  0 


(F- 


The  bracketed  term  m  equation  (F-8)  is  the  closed  loop  transfer 
matr.x  of  equation  (12),  and  consists  of  th*  numerator  polynomial  matrix 
1 Ni j (s 1 )divided  by  the  characteristic  polynomial  of  the  optimal  system. 


Jet 


,s.kNW  J.,e, 


Jett  Is-A  -  BHJJ 


However 


T 


N  (s)  H  e-e;  « 


which  yields 


det 


l)  1.  N  ,(s)h.  0  0 

1  ;lv  1  1? 


0  £.  N  . (s)  h  0  0 

l  .'i  '  i  / 

0  Ei  N?i(s)  hi>°  0 


o  r.  n  .($)  h.  o  o 

i  ui  i? 


kll.N  . 
i-  U 

(S)h.; 

A 

(S) 

k£.N„ . 
1  €  1 

(s)h. 

J  l? 

A 

(s) 

kE.N, . 
1  - 1 

(s)h1? 

A 

(S) 

kE.N.  . 

l  4 1 

c,)hi2 

A 

Cs) 

0  0 


0  0 


s  0 

0  s 


fi¬ 


lm 
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Appendix  G 

Open  Loop  Transfer  Matrix 


A  Definition 


M(s)  *  (is  -  A)  v  B 
B  A  i  rf  ram*1  Ma  trices 

-X.. 


X 


-  g  cos  v0 


A  - 

z 

u 

Zw 

8sin  '  Zw  V’ 

0 

0 

0 

0 

1 

Mu  *  Hi  = 

u 

•»M  •»  M  Z  i 

w  w  w 

•  8si„ 

Vo) 

M  • 
l| 

X6E 

X6T 

x«l» 

B  * 

-Z6B 

"Z6T 

"Z6F 

0 

0 

0 

M...  ♦  M 

6b  w 

Z6E 

M6T  +  Mw  Z6T 

M.c  4  M*  Z.E 
.  oF  w  6F 

Substituting  values 

from  Appendix  (A) 

for  the  F-8C 

-0.0506 

0.0072 

-33.8 

0 

A  = 

0.264 

-0-554 

127.4 

0 

0 

0 

0 

1.0 

-5.83  X  10 

-4 

0.00465 

-1.051 

-0-390 

-0-672 

0.001595 

-6.72 

1 

B  = 

19.95 

2.194  X  10'4 

18.9 

0 

0 

0 

-2.196 

6 

-4.67  X  10' 

0-0357 
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TRANSIHR  MATRIX  ELEMENTS 

Ma  (*)  *  C  S3  ♦  0 -730s2  -  109. 5s  -  53-9) 


Ml,.  («| 


Mi?  (s)  = 


Mu  (s)  ♦ 


Ml,  j  (s) 


-19.95 

aTTT 

(s?  ♦ 

2.20 

IT.:; 

fs2  ♦ 

2.20s 

Ms)” 

(s2  ♦ 

-0  0016 

,r  c  ^  4 

A(S) 

0.000219 ,  , 

M») 

— - 

4.67  X 

i£L6(s; 

MM 

4.67  X 

L°-6sfs 

Ms) 

s  1.S 

6.72 

Ms) 

Cs3  ♦  1 

-18.9 

Ms) 

Cs3  ♦  i 

-0.0357 

r<2  ♦ 

Ms) 

Is  ♦  , 

-0.0357 

Ms) 

s(s2  ♦ 

s4  ♦  0.995  s3  • 

0.021) 


18s  -  0-0218) 
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Appendix  H 
AMOAC  Design  Data 

System  0 

A.  Performance  Index  Parameters 


»  ^ 

b-  =  1.0  X  io'13 

b2  -  1.0  X 

10-13 

+7 

> 

4.10 

kj  =  4.0  X  10 

.  k  2  =  1.0  X 

10 

B.  Feedback  Gains 


Airspeed 

Sink  Kate 

Pitch 

Pitch  Rate 

Elevator 

>0.000282 

-0.000807 

-0.0679 

-0.1259 

Thrust 

1451.0 

■4977.0 

13,590.0 

222,600.0 

DLC  Flap 

0.0000199 

0.0000876 

0.00671 

0.00892 

C.  Numerator  Zeros 


N  (s)  =  kx  (s  ♦  1.408)  (s  ♦  2.352) 

W- 

u  (s) 

N  (s)  =  k2  (s  +  0.220)  {s  +  0.494  +  1.091j)2 

6TCs) 

h  (s) 

N  (-J)  =  k3  (s  ♦  2.86)  (s-‘.  12C  +  0.66j)2 

h  C*J 

N  (s)  =  k4  (s  +  1.68)  (s  +  0.119  +  1.403j)2 

6F(s) 

NADC-SD-7153 


Systea  I 

<\  Performance  Index  Parameters 


A 

1.0  x  10*13 

b, 

b? 

=  1.0  x 

I0*l? 

K 

4.0  X  10 

1.0  x 

*10 

10 

B  Feedback 

Gains 

Airspeed 

Sink  Rate 

Pitch 

Pitch  Rate 

U levator 

-0.000278 

— 

-0. 000794 

-0.0669 

-0.124 

Thrust 

1433.0 

4918.0 

9653  0 

214,200. 0 

OLC  Flap 

0.000196 

0.000866 

0.0663 

0.0877 

C.  Numerator  Zeros 


»(s) 

6b(s) 

N  (s)  =  kt  (s  *  1.411)  (S  +  2.324) 

U  Is) 
0T(s) 

N  (s)  =  k.»  (s  *  0.235)  (  S  <  0-493  ^  1.09j) 

<STCs) 

N  (s)  =  k3  (s  ♦  2.86)  (S-0.116  ♦  0.658)? 

h  (s) 
tfW 

N  (s)  ki,  (s  ♦  1.673)  (s  ♦  0-127  ♦  1.39812  j 

. .  i 
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System  3 

A,  Performance  Index  Parameters 

- - - , — 

b '  -  1.0  x  10"  13 

ki  -  4.0  x  10+ 

B.  Feedback  Gains 


1.0  x.  10. 


=  1.0  x  10 


Elevator 


Thrust 


Airspeed 

f  -0.000228 


1255.0 


DLC  Flap  0.00160 


Sink  Rate 
-0.000649  ‘ 


0.00756 


Pitch 

•0.0560 


0.574 


Pitch  Rate 


To  a.7" 


150,600.0 


u  (s) 

1  2 

N  (s)  =  k;>  (S  +  0.359)  (S  +  0.480  ♦  1.086j)  . 

•  i 

Last 

4>vCs) 

i  i 

N  (s)-  =  x3  (S  ♦  2.839)  (S  -  0.089  ^  o! 632j 

h  iil 

6p(s) 

i 

i  1  .  . . . 

N  (S)  =  kt*  (S  +  1.570)  (S  +  0.169  +  i.342j) 

*  J  »  1 
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System  4 

A.  Performance  Index  Parameters 
bi  =  1.0  x  10"l? 

k.  =  4.0  X  10+/ 

B.  Feedback  Gains 

Airspeed 


b2  *  2.0  x  10" 


k2  =  1.0  x  10 


Sink  Rate  Pitch  Pitch  Rate 


t levator 

-0.0001894 

-0.000538 

-0.0476 

-0.0866 

Thrust 

1116.0 

4209.0 

-17,030.0 

103.106.0 

DLC  Flap 

0, 00257 

0.0134 

1.007 

1.212 

C.  Numerator  Zeros 


ii£L 

Vs) 


«_(*! 

6"p  (s) 

Lt*i 

«T(») 


J*  J»). 

6  FCsJ 


N  (sj  =  ki  (S  ♦  1.6b4  ♦  0.210j) 


N  (s)  *  k2  (S  +  0.475)  (S  +  0.467  +  1.080j)? 


N  (s)  =  k3  (S  +  2.818)  (S  -  0.065  ♦  0.613j)‘ 


H  (s)  =  ki,  (S  ♦  1.478)  (S  ♦  0.201  +  1 . 294 i ) 2 


-4 
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System  6 

A.  Performance  Index  Parameters 

'  .  - 13 

bL  5.0  x  10 


k,  =  4.0  x  10 


B  Feedback  Gains 


b?  =  2.0  x  10‘ 

k?  =  1.0  x  10' 


Elevator 


Thrust 


Airspeed 


•0.0000465 


530.6 


Sink  Fate 


-0.000112 


Pitch 


-0.0134 


Pitch  Kate 


■0.0249 


DLC:  Flap  0.00691 


C.  Numerator  Zeros 


3070.0 


0.0018 


■2312.0  -54,380-0 


3.579 


3.000 


6,Csl 


6T(s) 


N  (s)  *  ki  (S  ♦  0.963)  fS  +  2.131) 


.  N  (s)  =  k;.  (S  +  1.600)  (S  +  0.336  +  1.063j) 


-  N  (s)  *  k3  (S  +  2.538)  (S  ♦  0.099  +  0.549))” 


~ N  (s)  =  kt,  (S  +  0.913)  (S  ♦  0.327  +  1.049j)? 
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Design  Limitations  of  Root  Square  Locu 


Technique  for  the  Carrier  Approach 


In  order  to  gain  satisfactory  phugoid  loci  in  the  outer  loop  closure 
of  the  root  square  locus  design,  it  is  necessary  to  have  a  lightly  damped 
complex  root  pair  in  the  sink  rate  loop.  This  complex  pair  is  obtained 
from  the  complex  zeros  of  the  sink  rate  to  flap  and  the  sink  rate  to  thrust 
transfer  functions.  A  method  of  testing  these  transfer  functions  for  suit¬ 
able  phugoid-control  zeros  is  derived  in  this  appendix. 

A.  Development  of  Phueoid-Control  Zeros 


For  the  derivation  of  N23(sJ  and  N22(s),  equation  (2)  is  simplified 
as  follows:  Since  ya  is  approximately  4  degrees,  cos  y0  is  assumed  equal 
to  unity.  It  is  further  assumed  that  a  suitable  mechanical  flap-to-ele- 
vator  interconnect  will  cancel  the  M^.  effect.  The  following  inequalities 
are  properties  of  the  aerodynamics: 

I  Z  |  *  •  g  sin  y0  (a) 


|  M  a  w  |Mw  2  a  >  ?  I  %  g  sin  Y  0  ! 


M  >  M. 
q  a 


Table  VII  lists  the  values  of  these  parameters  for  the  F-8C,  the 
A-5A,  and  theF-4t.  Inequalities  (a)  through  (c)  are  used  to  further 
simplify  equation  (2).  Thus,  the  airframe  equation  becomes: 


NADC-SD-7153 


*1  -z. 


-X  X  -g  0 
w  o  ® 


Z  -Z  Vo 
w  w 


0  0 


M  +M.Z  )  -M  +MwZ,  M 
u  w  u  w  w  a 


°  *  +  -Z61: 


;,6E+V  M6T+>\  MwZ6F 


♦M' Zxt.  H  +M*  z 
w  6L  /  ^  w  3T/ 


x  =  A  x  ♦  B  6_  (I- 

Laplace  transforming  equation  (1-1)  and  solving  for  x  /ields: 
/  1  / 

X(s)  =  [Is  -  A  ]'  B  A(s) 


T 

Airplane 

1 

z 

a 

g  sinyo 

M 

a 

M-  Z 
w  a 

M 

M* 

a 

F-8C 

125.2 

2.25 

1.074 

0.0237 

0.347 

0.043 

F-4B 

83.6 

2.25 

30 

•H 

0.0202 

0.433 

i 

0.0789  ; 

A-5A 

146.7 

2.25 

0.88 

0.0678 

0.288 

0.103.  I 

1 

F-111B 

iiiai 

mm 

0.08 

0.041 

0.436 

0.0856  j 
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Ci2fs)  -  -Z  (s2  -Ms-  Mo)  -e  z  Vo 
1  *  •  '  u  q  w 


The  magnitude  of  Z  Mot  and  S  Z  Ve  are  compared  for  the  selected 

U  w 

aircraft  in  Table  VIIJ. 


Airplane 

8  Zw  Vo 

Z  M 
u  a 

F-8C 

0-073 

0.284 

F-4B 

0.032 

0.396 

A-5A 

0.00683 

0.0479 

F-111B 

0.00743 

I  0.0247 

TABLE  VIII  CONSTANT  TERMS  OF  C12 (s) 

This  comparison  a3 lows  the  simplification 

C 12 Cs)  =  -  Zu  (s2  -  M  s  -  Mu) 

H 


(I-5a) 


s  -  Xu 


C22(s)  =  det  0  s  -1 

-0  -Ma  s  -  Mq 

Cj!;>(s)  =  s3  -  (Xu  +  Mq)  s2  ♦  (Xu  Mq  ♦  Ma)  s  +  Xu  Mj  +  (g~xa)  e 

Tne  magnitude  of  terms  in  the  coefficients  are  compared  in 

Table  IX.. 


Airplane 

Xu 

Mq 

XuMq 

it  i'  • 

M 

« 

XuMo 

(8-VB 

.  F-8C 

0.0506 

0.347 

0.0176 

1.074 

0.0543 

0.0195 

i 

F-4B 

0.094 

0.433 

0.0407 

1 . 298 

0.122 

0.00759 

A-SA 

0.0544 

0.288 

0.0157 

0.83 

0.0479 

0,000914 

F-I11B 

0.0456 

0.436 

0.0199 

0.08 

0.00365 

0.00117- 

TABLE  IX  COEFFICIENT  TERMS  OF  C72(s) 


'-'-VM 


ifTW 


NADC-SD-7153 

Thu*  C  (s)  is  Approximately 


C22(s)  -  s3  - 

Mq  s2  -  Mas  ♦  Xu  Ma 

(I-Sb) 

s  -  Xu  g  -  Xa  0 

C42  (s)  =  det 

Zu  Zw  Vo  0 

0  s  -1 

642(3)  “  -zw 

Vo  s  ♦  ZwVoXu  +  zu  (g  -  X#) 

•  ] 

(I-Sc) 

a  Sink  Rate  to  Thrust  Numerator  N22(s) 

4»  '  > 

4 

N22(s)  *  2  C^2(s)  b^2 
i=l 

N22(s)  b  -Z6t  S3  ♦  (Z6lWq  -  ZuX6x)  s2+ 

[  2U  Z6t  Mq  ♦  ZfiT  (Ma  -  ZwVoM;)  -  ZwV„M6t]  s  ♦ 

(Z„  MaXjT  -  z{T  Xu  Ha*  C*-X„)  zu  CMST  *  m;  Z„)  * 

*  W.  <*{T  *  K  W  l 

Table  x  lists  the  magnitude  of  the  various  coefficient  terms. 
N22(s)  may  be  simplified  to 

N22(s)  =>  -  ZfiT  [s3  +  Z^s2  ♦  (ZwVoEt  -  ZuMqr,T  -  M0)  s  -  Ma  nT  Zj 

.  (1-6) 

for  all  aircraft  listed  in  Table  x  except  the  F-111B,  nT  and  eT  are 
defined  as  follows: 

X*T  M-  ’  ■ 

and  et- 

The  extremely  small  value  of  M0  for  the  F-111B  negates  the  constant  term 
approximation  of  equation  (1-6).  This  approximation  will  be  retained  and 
justified  later  on  for  the  F-111B. 
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For  the  F-8C,  N22(s)  has  a  complex  root  pair  in  the  right-half  s 
plane.  Inspection  of  N?200  in  Table  xi  shows  that  the  quadratic  tern: 
formed  by  deleting  the  s*  term  also  yields  a  complex  root  pair  in  the 
right-half  s  plane.  Table  XII  list*  the  roots  of  N22  (s)  when  this  quadratic 
term  is  assumed  as  a  factor  of  N?.?.  (s)  and  also  when  the  original  cubic 
equation  is  solved. 


! 

Airplane 

Quadrat ic 
Approximation 

Original 

Cubic 

F-8C 

♦  0.25  +  1 -01 j 
-2.39 

+  0.34  +  0.82j 
-2.59 

A-5A 

-0-19  ♦  0.92j 
-12.20 

-0.16  +  0.96j 
-12.28 

F-4B 

+0.07  +  1 . 14 j 
-3.43 

+0.20  +  1.07j 
-3.74 

F-111B 

-1.71  -  .05 

-2.30  -  .047 

-7.35 

-6.76 

TABLli  XII  ROUTS  OF  N?2  (s) 

Using  the  approximate  quadratic  factor,  the  phugoid-control  zeros 


of  N)<(s)  are  described  by 
'  /  ZwVoM 

s  ♦ 


oT 


MaZ6T 


Zu  X6T 


-  Mq 


) 


M 


(1-7) 


For  the  phugoid  locus  to  pass  through  the  desired  region  of  figure 
12,  the  numerator  roots  should  have  a  natural  frequency  greater  than 
0.707  radians  per  second  and  a  damping  ratio  no  greater  than  j  0 . 35  | . 

Thus,  the  phugoid-control  zeros  which  determine  the  numerator  roots  should 
have  a  damping  ration  between  -0.35  and  0.35  and  a  natural  frequency 
greater  than  0.707  radians  per  second,  resulting  in  the  following 
specifications: 
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Table  XIII  lists  the  magnitudes  of  the  coefficient  terms.  Dropping 
the  less  significant  terms  yields 

-  Z&¥  [S3  ♦  (Zu  np  -  MqJs?  -  MfS  -  Mn  nF  -  X0)]  (I-9J 

where 


Table  XIV  lists  N2j(s)  for  the  F-8C,  A-5A  and  F-111B.  The  F-4B  is 
excluded  due  to  the  nonlinearity  of  its  flap  effectiveness  cuxve. 

Inspection  of  N23in  Table  XIV  discloses  a  difference  of  two  orders 
of  magnitude  between  the  linear  term  in  s  and  the  constant  term.  This 
fact  strongly  indicates  that  it  may  be  approximately  factored  as  follows: 
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N,*(s1  r  -Z.c  ts-'  *  (Z  nF'-  M  -)s?  -,M  s  -  M  (Z  nF  X  ) 
6F  *•  .  u  q  1  >  a'-u  :  u-' 
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lhus  the  phugoid  control  zeros  of  N; j  (s)  may  be  described  by 

,  *  (Zy  -  M^)s  -  Nb  =  0  (1-10) 

bo t  the  roots  of  equation  (1-10)  to  have  a  damping  ratio  between  -0.35 
and  0.35  and  a  natural  frequency  grearer  than  0.5  radians  per  second,  the 
following  inequalities  must  be  satisfied: 


-M  (V 
a  — 


2  X  -  M  Z... 
U  SL-IL. 


f 1 • < lb ' 


P.  Operational  limitations 

To  assure  closed  loop  phugoid  roots  which  satisfy  equation  (2o.  a 
lightly  damped  complex  zero  pair  is  needed  in  the  sink  rate  loop,  ^in  » 
tin  elevator  zeros  are  characteristically  real,  this  pair  must  come  tV'ir 
either  the  flap  zeros  or  the  thrust  zeros.  Alternatively,  the  airodynamic 
derivatives  of  the  airplane  in  question  must  satisfy  either  conditions 
(1-8)  or  conditions  (.1-11).  Since  conditions  (l-8a)  and  (I-lla)  are 
identical,  it  appears  that  -M^Q.S  is  an  essential  requirement  for  satis¬ 
factory  closed  loop  phugoid  roots.  Table  XVF  lists  the  results  obtained 
when  these  conditions  are  applied  to  the  F-8C,  the  F-4B,  the  A-5A  and  the 
F-111B,  each  for  a  single  approach  configuration.  The  results  indicate 
that  the  Hoot  Square  Locus  technique  can  be  applied  to  design  satisfactory 
automatic  approach  controllers  for  the  F-8C,  F-4B,  A-5A  and  similar  air¬ 
craft.  In  the  case  of  the  F-111B,  however,  the  very  low  magnitude  of 
Ha  indicates  that  adequate  phugoid  bandwidth  cannot  be  obtained  via  root 
square  locus  design 
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This  particular  value  of  -M  -  008  is  for  a  specific  loading  condition 
as  shown  in  Table  XV r.  Examination  of  the  entire  range  of  center  of  grav¬ 
ity  locsTi.'pc  and  airspeeds  of  the  F-111B  approach  configurations  indicate 
that  -v!ue  of  Ma  is  always  unsatisfactorily  small  This  may  be  dem- 
oii'-tnted  as  follows: 


S  cp  ; 

V0?  C  aid 

2ly  tv 


C  =  -  C 
m 


/xn  xc  ) 

L  (“I - ^3 


Substituting  in*o  (T-13) 

/s  c.  \ 

\  \lTy — )  V"  ' 


(X  -  x  ) 


(1-13) 


(I-’.’.i 


(1-14) 


Airplane 

Landing 

Conditions 

F  8C 

Vc 

=  22 6  ft /sec 

tg 

*  27%  MAC 

F-JB 

Vo 

-  22b  ft/sec 

*  30-0%  MAC 

A-5Z 

Vo 

=  223  ft/sec 

‘-'g 

=  30.6%  MAC 

F-111B 

Vo 

--  200  ft /sec 

‘-'g 

-  35.5%  MAC 

1  074 
1 . 298 

oTslT" 
0  08 


A 

o.4y4 

0.21O 

1)7378" 

i  7:77" 


u 

0.25? 

Not  cal¬ 
culated 

0  231 
"77374 


NOTH ;  Satisfactory  zeros  require 

-M  •  0.5,  A‘0.5  and  0  0-7 

a—  —  — 

TABLE  XVI  CRITERIA  FOR  SATISFACTORY  PHUG0 ID-CONTROL  ZEROS 

In  the  approach  configuration,  is  independent  of  the  center  of 
gravity  location  and  airspeed.  Therefore 


Ma  -  k.  Vo^  (x„  -  x„  ) 


(1-15) 
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ar.J 

Vo  ■•''50  ft; sec 
which  yield 


M  •  1  ’b  * 


|M0|  *  0-332 

Thus,  the  magnitude  of  Mq  for  the  F-1UB  is  unsatisfactory  for  any  center 
of  gravity  location  and  approach  airspeed,  in  terms  of  root  square  optimal 


design. 

Condition  (l-8a)  was  derived  from  the  approximation  of  defined 

in  equation  (.I-bJ .  This  approximation  was  f6und  to  be  poor  for  the  F-'lllJ 
for  the  flight  condition  given  in  Appendix  A.  Several  components  of  the 
constant  term  are  similar  in  magnitude  to  n.p  Zu  Mft  at  this  flight  condition. 
However,  for  the  power  approach  conditions  stated  in  equations  (1-17 >  a.’! 


of  the  components  of  the  constant  term  except  n,^  Zy  Ma  remain  near  the 
values  listed  in  Table  x,  while  n.^.  Zy  Ma  increases  by  a  factor  of  five. 
Thus,  a  natural  frequency  of  sufficient  magnitude  (i.e.,  =  0*07)  can 

only  be  provided  by  the  Zu  M0  term.  Conditions  (I-8aj  is  thus  a  valid 
criterion  for  the  F-111B  airplane. 
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